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ABSTRACT 


The objective of this thesis is to investigate the problem of the 
Safe service J/fe prediction of graphite composite structures such as 
аг га меззейз usec in tne pilot eject!on seats, rocKel motor cases and 
space shuttle energy storage compartments. 

The basic data required for life prediction is the stress rupture life, 
i.e., the composite life under constant load and it has to be recognized that 
extensive testing over /ong periods of time is required (0 ргоачсе 
statistically meaningful data. The contribution of this investigation 15 
focused on identifying the characterization methodology for efficient data 
generation (to minimize cost and time), the appropriate theoretical models 
for correlation of the data (to translate data to applications) and in 
evaluating tne limitations of the available testing (lower bound response 
speed) and data acquisition equipment (upper bound recording rate) thereby 


enhancing future experimental planning and testing. 
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I. INTRODUCTION 


А. BACKGROUND 

Although composites have been employed ever since materials were 
first used, such as the naturally occuring wood or the man-made mud bricks 
reinforced by straws and have been widely used in the recent 50 years In the 
building industry (reinforced concrete), the high technology of tne 
composites has evolved in the aerospace industry only in the last 25 vears 
The first strength critical application fer modern composites were 
filarnent-wound pressure vessels using glass fibers. 

What has highly motivated the effort of composites development 
has been the goal to come up with a combination of properties not 
achievable Dy any of the constituent materials acting alone. Thus a solid 
could be fabricated from elemental materials which oy themselves could 
ши лети претисшаг Gesign requirements Such as strength and/or 
eee ss lO Volga га оза se Of WHICH їз OT Crustal importance in the 
aerospace industry. Another highly motivating aspect was the very їр 


reliability of the cornposites due to the micro-redundancy which.cornes as 3 


result of the fiber-matrix load sharing that in fact diffuses the local fiber 


failure sites. And it was this materials redundancy that finally proved that 


composites were a superior application for monocoque structures such as 


pressure vessels and rocket motor cases. 


According to Tsai [ Ref. | ] the results of this effort in the fielen 


the Aerospace Industry came sequentially in the form of vartous fuselage 


components as well as sorne primary control surfaces of several first class 


Hanter: ЭЦ: 


Т. 


The boron/epoxy fuselage section and horizontal сат ора 
General Dynamics F-1 М. 


The graphite/epoxy fuselage components for the Northrop F-5 
made also by General Dynamics. 


The limited use of the boron/epoxy material system in the rudder 
of the Mc Donnell-Douglas F-4 which has accounted for a 558 
weight reduction compared to its Alluminum counterpart. 
The boron/epoxy horizontal stabilizer of the Grumann F- 14. 


The composite stabilator of the Мс Donnell-Douglas F-15. 


[he graphite/epoxy horizontal and vertical stabilizers of the 
General Dynamics F-16. 


On the other hand in the commercial aerospace industry Boeing was 


the first manufacturer to utilize about 2 tons of composite materials in the 


10 


floor beams and all the contro! surfaces of the 767 airliner while the Beach 
Aircraft s Starship | all-composite aircraft is scheduled to start its 
flight test program before the end of 1996. 

Apart from the U.S manufacturers extensive use of composites has 
also been used in the Dassault-Breguets Rafale, the Israel Aircraft 
Industries Lavi and the all-composite fin box of the Airbus Industrie 
небо an impressive-structure in its simplicity, also described by Tsai 
[Ref. 1] 

In paralle! to the manufacturing milestones that had to be resolved 
once the goals for composite compoments were set, the reliability proolems 
had to be resolved as well. Among them was the problem of safe life 
prediction of the composite components and its relation with the composite 
Strength. The significance of solving this problem for the aerospace 
Industry becomes more obvious when one focuses on the construction of 
Space and/or missile applications where repetition for experimental 
purposes is either impossible or limited to numbers often less than 3 due 
to obvious economic constraints. Consequently one has to realize that other 
methods have to be established, methods such as statistical inference and 


otner mathematical techniques which will enable the engineers to estimate 


accurately the service life of the composite component or system under 


TH 
СУ 
iD 


evaluation, based on a limited number of experiments which I3 accep 
from the cost and time standpoints. 

From physical experiences and mathematical modeling, It 1s Known 
that the average macroscopic strength of composites depends strongly o 
the number of micro-fiber failure sites and the strength of specific 
composites depends on the clustering of the fiber failure sites. AS amom 
information on the statistical characteristics of fiber filament strength 15 
the base for estimating the reliability of the macro-strength composite 
structures. Specifically the statistical! tower tai- distribution спее 
filament strength and filament life are nad The characterization of the 
lower tail behaviour requires testing of large number (thousands) of 
samples. In the case of strength characterization, large number of samples 
Is.a matter of economics whereas in the case life characterization, WENNS 
matter of economics and time. Aside the cost constraints which are always 
obvious, the reason that time constraints are rnore than equally seriouamin 
, 18 that experimental testing of fiber filaments could last 
ingefinitely without providing adequate data. This inherent uncertainty in 


lire characterization can only be overcorne by acceler — testing 


29 


methodologies. In this investigation, fiber bundle testing is chosen to 
increase the number of samples. However, fiber filarnent testing IS si! 
required to serve as a validating tool on the fiber bundle testing results. 
Higher stress levels are used to accelerate the life testing Іп order to 
assure that failure data will be observed within the existing time 
constraints. Under stress-life acceleration circumstances the service life 
prediction could be made possible at a given confidence level utilizing tne 
appropriate statistical methods, given that the fiber bundles under testing 
nad already been able to withstand either a prescribed load level without 
any failures, or a higher load level experiencing a Known number of filament 
failures per bundle. The results of this investigation are relevant to 
man-safe applications as in composite pressure vessels used in the pilot's 


Sueetion seats, the liquid oxygen tank in the jet aircrafts, the space shuttl 


ср 


tanks for energy Storage, or the cornposite rocket motor cradel of the 
guidance system of various missiles, where mostly graphite fibers are 
used, as well as the broad class of composite applications where the 
SupsicLume. i5 subected to sustaimed tensile-deading. 

Пре prooler of accurately determining the life of one of current hign 


strength graphite fiber (Hercules A54) has already been started at the N.P.5 


21 


by Lt.Fred D.Carozzo and tme results Oj маг s[Ref. 2] identified three 
[rue e PNIS SRI TIED that have to be = игајпајед in order го mars inis 


characterization methodology operational: 


й To avoid the stalling of the tensile testing equipment during tne 
loading sequence, because Stalling appears almost exacliy 25 ап 
increase in broken fibers in a bundle and hence can produce false 
interpetation. 


2. To explain why the creep dala recorded during the creep-rupture 
phase of the experiment do not agree with the expected phvsics 01 
the problem. 


E. To increase the recording sace of the data acquisition system 
used, because it is a known fact that much more fibers per Dun 
must have failed during the loading sequence than what has 
already been recorded. The number of fibers per oundle within tne 
bundle which failed during ine loading sequence Is needed 10 
calulate the conditional probablity for the A54 fiber bundle 
servi Gems. 


EI 
ср 


D. GOALS 


The objective of this researc en (ст? 5? от өт. 


~ 
A 

i 
9 
т 
23 
to 


reliability characterization methodology is to obtain strengtn- HT eco NEN 
grapni*c fibers. The final target is to present these data in Sücn a мани 


"Given a8 stress level one could predict the probability of Заг сре 


~ 
ha 


for a known confidence level and vise versa (that is, given the anticipated 
service life of a composite application specify the stress level that snould 
not be exceeded during service in order to achieve the given safe-life limit 
within a prescribed confidence level) ". With this ultimate goal in mind, 


what this thesis is desired to contribute to is four fold: 


|. To present a general lay out of the methodology that is to be used 
as a guideline to reach the ultimate goal of obtaining strength-life 
data. 


2. To evaluate what is the minimum loading speed that the existing 
testing equipment can be operated with, so that no stall will occur 
during the loading sequence of the fiber Dundles in order to avoid 
interruptions and to identify the load dependency of the .minimum 
speed. 


3. To give a reasonable explanation on why the creep data recorded 
during the creep-rupture phase of the experiment do not agree 
о г:20 025105 ОГ the problem. 


4 To study the limitations of the existing data acquisition system 


and propose solutions to increase the recording rate of data points 
during the short time loading sequence. 


Consequently the projected end-use for the information generated bv 
this investigation can be specifically tailored to the applications of the 


Hercules AS4 graphite fiber bundles and hence. the upper load level 
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boundaries for specific applications can therefore be established such that 





the predetermined service liie ef the composite component under 


evaluation could be safely achieved. 
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|. APPROACH TO THE PROBLEM 


A. GENERAL 

It is intuitively obvious that due to the minute transverse dimensions 
of the composites fibers, huge numbers of fiber bundles have to be used in 
ШЕ 051206 10701! а согрозме Structure, regardless of its absolute 
dimensions. It 1s also evident that the type of constituents to be used in a 
composite structure ought to have been experimentally tested in advance 
and demonstrated that they comply with the desired standards and/or 
ввеснесасіов5 сопсептіпд their strenqtn and service life prior to the 
beginning of the manufacturing process. This means that there has to exist a 
pilot model design fabricated under Known parameters prior to the 
ШШШПЕПСШЕҮЕП: ОСІ tae  meanublaciuping process. Furthermore, since the 
constituent materials used for the model evidently differ in essence from 
those the actual manufacturing will start from (even though they have been 
fabricated from the same material compositions and according to the same 
procedures) there is a need to assure the reliability of the production line 


components. The answer to production reliability can only be given in terms 
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of confidence levels and under agreed upon conditions between the Procuring 
Activity and the Manufacturer. It is also likely to assurne that within the 
oroduction the fiber bundles will not have the same strength or the sarne 
life. Heuristicallly, it is assumed that " //gh strength is associated 
with long life whereas low strength with short 
life respectively ^". Under the above circumstances, provided it is 
possible to identify the strongest and the weakest fiber bundles by quality 
control methods, Figure ! could serve as a good starting point for 
understanding how strength and life of the same fibers could be related, 
under the proviso that these relations have to be determined as accurately 


as possible later on. Rosen [ Ref. 3] and Phoenix & Wu [ Ref. 4 ] show that 
both the strength (at t = tọ where Р, is a very short finite period of time), 
and the life of the fiber bundles at various levels of strength follow the 
Weibull distribution which is usually described by means of a snape 


parameter (ar) and a scale parameter (9). 
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Figure 1. ‘Weibull Strength to Weibull Life correspondence 
of strong and weak fiber bundles composed of 
iQ fiber filaments. 
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В. BASIC LINE OF METHODOLOGY 

Without any loss of generality one can now follow the  inouüan: 
experiment described below for a fiber bundle composed of ten (iO) 
filament fibers identified from now on as f, for i=1,2,3,...,10. One can also 
assign the intrinsic strengths of these Filament fibers to be (with «Пе 
implicit assumption that the numerical values of tne strength random 


variable follow the Метли distribution) as follows: 


ii 


1. One filament (assume T of SEED РА СИН КОЈИ А 


гэ 


i 


Опе filament (assume f5) of strength Pe - 12 108. 


ій ( - | Г ( - 
3. Two filaments (assume f^ and f) of strength В. 
4. Three filaments (assume fs, fa and i7) of strength Pc 4 =15 155. 


5. Two filaments (assume fg and Га) of strength Pas = 18 Ibs. 


о. One filament (assume f 49) of strength Pag 5 20 105. 


f one desires to estimate the Strength of this bundle, one Сайн аа 
thc bundle to different strength levels ( P1, P2, P5, etc.) and оре си 
survival of the ftlaments. There wil! 5e strength levels belosve ПОП 


~ 


the weakest fiber of the ten in the bundle can safely withstand. For a 
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continuous increase in stress level there exists be a load level Ша ја 
which the weakest fiber will Tinaliy jail and, under the same (лоцогі 
process, there must be another load level at which the next weaker Tiber in 
the bundle will fail 22) and so on, until all the fibers in the bundle wii: 
fail and hence one will be able to observe the failure of the entire buncie., 
Having already postulated that the fioer’s strength in the sundie faraws 


the Weibull rnodel distribution, one can plot the probability of failure for 


гр 


this bundle at t з Зара з presented in Figure 2. At this point it nas fo c: 


noted that the finite time tg is only of the order of a few minutes and 


Wwemerore- negligible if compared with the anticipated duration of сле 
service lite and therefore, with tne experiment to foliow Гог the пе 
testing. 

A very important observation =o be noted is thal the fibers in ine 
шэг аге aucomatically ranked with respect to their strenain, which 
Пеле tie Weorer 1iders In the Oundle fall first whereas the stronger 
E 212 ог ас а higher load level. The significance of this 
observation is what actually explains the big advantage of the fiber Oundle 


testing versus the single filament fioer testing and can become more 
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Figure 2 Homologous correspondence of мери! 5 гепатл а Е + 10 
to Weibull Life distribution for afiber bundle composed of 
10 fiber filaments. 


obvious if one thinks of bundles composed of 5000 up to 10000 filaments 
(as it is usually the case). the expediency is (лас ше lower tali wii 


always be the first to бе observed. As it has already been stated earlier, 
if the applied load level does not exceed for example, the (Py |) leve], then 
no failures will occur until n However if one were to leave the bungie 


«р 


У) 
w^ 


гр 
шор 
с 
т 
с) 
гэ 


жа ныға ғ. и : ^ a e 4 
some failures couid be observec 


с 


long period of time. And it is reasonable to think that under tnis load Р. 


after a long period of time the weaker fiber will fail first, the next weaker 
fibers will fail later, and so on, until all the fibers in the bundle will nave 
failed, when it would also be possible for someone to observe the failure of 
the entire bundle. Assurning again a Weibull distribution model for the lire 
of the fibers bundle as per Phoenix & Wu [ Ref. 4 ] one can plot the 


ВА oF failure of this very buncie in time or in other words tn life, 


given that no failures have occured during the loading sequence until (o 7 


я ща с - қа В чини ыша 


* How this plot could be obtained will be explained later in ine 
Detailed Method of Solution. 


Usually due to very long period oí Lime which may have to elapse until 
the first failures could be observed it 15 озса о ЗЕ 
Log(t) coordinates as far as the abscissa is concerned апа из 18 
the way the figures will be presented during this thought experiment. 

But the situation described so far is merely the simplest, and 
therefore the thought experiment nas to be carried through In a rnore 


realistic way. Such a real world application could 5e one in whicn tne ioad 


level (P,,,4) and this load level could nave been maintained to be constant in 


ум 2 


time corresponding to a stress level say (Qa). 


Under these circumstances what should happen could be described as 


follows: 
1 During the loading sequence until to the fibers E OMS > 2,1) агд 
fo (Po > Poo ) will fail, whereas the remaining (Гэ “үс 
Survive. 
2. During the time that the constant load P2 is maintained, the 


rest of the failures have tu 3ccur in such a way that the ДЕ 
fibers f= and f, will fall earlier, to be followed Dy the failures o 


filament fibers fe, fc and 12, the гапе Зао анале ти - 


=, [6 8 


the failure of the last one Pip which completes tne Ізішге бі 


entire bundle 


In this simplest case where no failures had occured during the 
loading ppm: (due to the low magnitude of the applied load) one was able 
to infer about the life of the fiber bundle in terms of the Weibull model 
based on the magnitude of the load level the fiber bundle was able to 
withstand without failure. Direct statistical inference was possible 
Æ Use no failures nad — during the loading sequence. However in the 
real world applications there are always finite number of failures during 
thc loading sequence and therefore one has to be able to determine tne 
Eu smolthese failures. and their respective failure stresses in order to 
be able to use the same technique and infer about life. In other words life 
e ence 15 possible only om the basis of the conditional probability 
HOC failures did not occur during the loading sequence. And in order to 
evaluate this condition one has to Know the exact number of the fiber 
filaments that have failed during the loading sequence, so that the 
magnitude of the load level that would have been achieved if no failures had 
occured can be used as before. This requires an accurate data acquisition 


Моа паз to De used in order to be capable to record tne possible 


failures during the loading sequence. 
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Although tne sequence of failures desenibed зо паг Бис 


anticipated in accordance with intuition, (i.e whatever strength of tne fiber 
filaments was not consumed in testing during the loading sequence, can De 
used for longer life endurance under a homologous strength-life 


relationship), this might not be the case and therefore one has to find out 


experimentally whether -this one to one correspondence exists between 


failures at t» t, and failures in time ( t - t; for 1 7 ab,¢,d ) Того ы 


0 
arbitrary selected load level or in other words whether the above 
mentioned relationship of strength and life is homologous or stochastic 
(see Figure 2). 

Оп the other hand, another major constraint is the fact that for all 
practical purposes it is not possible nor necessary to wait until the failure 
of the last filament fiber because (his might occur much later than the 
structural lifetime. This is the reason why the engineers only have to wait 
until tne first few failures occur and then terminate the experiments on 
purpose. The background of such a descision lies on the fact that the 
engineering Interest is іп the early failures of the structure 90088 


evaluation whicn justifies a safe service life. In terms of the thougnt 


experiment described so far if one assumes that the Ру? load level would 
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loading sequence at t = t,, and hence Tiber rilaments 


ad leve 


failed and furthermore the loa 


бе achieved during the 


Г, and fo would have 
managed to be reached in time, one should nave to wait at least until the 


first few failures (namely those of fiber filaments fz and fg) will occur 


when the experiment can be terminated. If the same thought experiment 1s 
reaching higher ioad leveis in everv ioadinj 


times 


carried out several 
is plotted every time as a Weibull distribution 


sequence and infered lite 
based on the load level that would nave been achieved if no failures nad 


happened during the loading sequence, a curve resembling the familhar S-N 


НЭЭ 
A 


fatigue curve will arise as shown in Figure 3. 
The significance of coming up with a graph such as that of Figure : 


for a real word application as the widely used araphite fiber A5-4 1s more 
than obvious, apart from the fact that serves as a basis for the so called 


Preof-/est with the intent of determining definitively beth the strength 


and the lite of the same specimen simoultaneously. 


Strength (PJ 


Figure 3. 





Life ( Log t? 


Weibull Life distributions for four identical 
fiber bundles comoosed of 10 fiber filaments 
acnieved as a resuit of different load levels 
reacned at the end of the loading sequences. 
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C. DETAILED METHOD OF SOLUTION 
What remained to be clarified is the method of representing the 
probability of failure of the fiber bundle in time, given that no failures nave 


occured during the loading sequence (O «t cta) where tẹ is a short period of 


time of the order of several minutes. 

To demonstrate how this can be done and also get more insight in the 
process of approach to the overall problem we will continue with the 
thought experiment considering four (4 ) identical fiber bundles like the one 


which has already been described. If one starts loading them up to four (4) 


different load 1еуе!5 namely Po = 101105, Ро2 3 12 165, Рез 7 14 Ibs and 


Ред” 16 lbs and then maintain these four load levels in time the following 
pictures can then be presented for each one of these fiber bundles: 


Ж CASE A. \0 filament bundle under P.,=10 Ibs tensile load 


61219 І05 (805116 1080 





since P.; 1s the strength of the weakest fiber filament in the 
bundle, obviously no failures will occur during the loading sequence. 
However, atter some period of time failures will start occuring as time 


increases, giving the picture shown in Figure 4-1a. Observing carefully the 
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steps of this curve it will not be difficult for one to understand that each 
ШИ 775: 2 2505 Stands fer certeir failures of fiber filaments in khe 
bundle. And the way this stepped curve is plotted shows that the fiber 
filaments of the same strength fall simoultaneously. S5moothening the 
stepped curve obtained, it is not difficult for one to observe that the new 


smoothened curve stands for the Reliability Cummulative Distribution 
Function (RCDF) of tne fiber bundle under evaluation. And since the Failure 


Cummulative Distribution Function (FCDF) is always related to the (RCDF) 


by the relationship: 
(11.1) (FCDF) = 1-(RCDF) 


it is evident that FCDF for this bundle can be obtained by plotting the 
complement image of the RCDF with T Lo the vertical (load ) axis as 
не 2-15. | has to be noted at his point that this picture is 
usually presented in the so called Weibull coordinates due to severa! 
reasons sucn as Che follow ing: 


a. The Weibull rnodel is used to 
Е 


со simulate both the strength and the 
life of the fiber bundles ([ Ref. 31a 


nd [ Ref. 4 |). 


2 
> 
2) 


D. When Weibull coordinates are used for life the abscisa 1s 
presented in Log (t) which is consistent with what has been 


already stated on the same aspect. 


ite 
(б 


3 
„1 
“4 


CD 


о 


In Weibull coordinates the FCDF represents itself as а Straignt 
line instead of a sinusoidal shaped curve as when physical 
coordinates are used. 


Гэ 


Therefore it is customary to present the picture of Figure «ёс! 
in the Weibull coordinates forrn as shown in Figure 4- 1c. Taking now into 
account that no failures have occured during the loading sequence, one can 
now plot the failures’ probability density function (fpdf) using the 


relalionsnio: 


(11.2) (Граг) - | d(FCDF)/dt] 


This failures’ probability density function for this first fiber bundle denotec 


2 CASE B: 10 filament bundle under Ре 212 lbs tensile do 


ые аны 


NJ 


Since Pe.» is the strength level that is higher trom onmyeene 


Nbe filament strength, it is reasonabie one to anticipate one laliure то 


occur during the loading sequence and if this load level 1s maintained in 
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time some more failures will start occuring until tre whole fiber bundle 
IE iepicture as the one shown in Figure 4-2a. The oniy 
difference between this case and the previous one is that the stepped 
curve obtained this time stands for the Reliability Cummulative 


Distribution Function of the fiber filaments that suryvsved the load level 
achieved at the end of the loading sequence (RCDF).. In order for one to be 


E cuseheésame technique as in the previous case (0 come up with the 
га шгез probability density funclion it Is first necessary Lo evaluate the 
load level that would have been achieved during the loading sequence under 
the Condition that no failures had occured. This load level can be easily 
avaluated given that tne number of the fiber filaments that have failed 


during the load sequence is known. If one then denotes this load level by 


PNF it is clear that: 


(1.5) (Реј /Рај ) = ( 100/% of survived filaments ) 


for I= 1,2 5 Ие 


Hence for the case under evaluation ( 1 7 2), where P.4- 1215s and the 
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percentage of the survived fiber filaments is 0.90 ( one failure out of ten ) 
one can easily find that EIU = 13.3 Ibs. Then the upper iert portion oi 
the existing stepped curve can be extrapolated to the left to account for the 
(RCDF)ye that is the Reliability Cummulative Distribution Function under 
the condition that no failures had occured during the loading sequence. Using 
the sare reasoning as before one cai now obtain the FCDF for both ins 
survived fiber filaments in this case (FCDF), and for the fiber bunble under 
the condition that no failures had occured during the loading sequence 
(РЕСОР yr. Of Course what is actualiy important for this case Is only ine 


latter that can De evaluated by equation (I.1). This result is presented in 
Figure 4-2b in physical and in Figure 4-2c in Weibull coordinates. 
Finally one can plot the fpdf for this case using equation (11.2) as presented 


in Figure 4-2d. 
D CASE C: IO filament bundle under edes lbs tensile load 
This case is exactly similar to what has already Deen described 
т О Пе 240 or 502552761 that would Nave deen 


acnieved at the end of the loading w.zquence if no failures had occured 


(pe... = 17.5 lbs) Therefore the whole procedure is outlined only 
schematically as shown in Figures 4-3a , 4-3b , 4-3c and 4- d 


4 CASE D: 10 filament bundle under Р, am 16 105 tensile load 


Again due to the fact that no conceptual differences existi no 


matter what is the creep-rupture load level the procedure is outlined oniy 


12 


schematically < У = 26.7 195), 23 shown in Figures 4-4a, 4-45, 


4-Ac and 4- 4d. 
Э. 22254//ё$ 
If one now takes the results of Figures 4-1d, 4-24, 4-34 апа 
4-4d and plots them with respect to tne corresponding load levels achieved 


at the end of each loading sequence, the picture shown in Figure 5 can be 


obtained in which the Sirength-Life curves for the extreme lower tal! 
(5-L),4 , the median (5-L)q,, and the extreme upper tail (5-L),, of the 


Weibull model distribution are plotted A close observation of (hit 1ig@ure 

can lead to the following worthy of corarnent thougnts: 
a. Given thal one can produce accurately several points on Lnemesus 
urves the exact determination of the life 415 04 ом 


z 
possible at any desired strcnath (stress) level is ӨЛІМ атпаса 
the availability of equiprnent if tirme is not a problem. 
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cose s Senec = LCife=curves for the extreme lower tail (S-t)It 
median (S-Ljm and upper tail (S-L)ut, based on the results 
‚ОГ four exoeriments of identical fiber bundles composed of 
lO fiber filaments under different Creeo-2upture loads. 


47 


The higher tne number of fiber filaments in the bundles The 

піапег the confidence level of the resulis го 08 обгатей, ог и 

otnerwords the more certain the ХУ? curves will De ani 
| 


The Weibull curves are obtained on the basis that one can produce 
an accurate curve FCDF under the condition that no failures nag 
occured during the loading sequence. The accuracy of the plot 
mainly depends on the exact evaluation of tne fiber filaments that 
had failed during the loading sequence and on the assumption that 
every fiber is carrying tho same amount of 'cad, 2c сана 
internal friction is absent known as the "egua! /oad-sharing 
rule 


іл 


What is of great practical importance is the exact determination 
of the extreme lowers-taH.o: rne Weibuli diztrioUclons Гог гаш ш 

load (stress) levels. Becallise this Ts onen 271504 61.4 -—— 
the safe service life of any application. Another very serious point 
in discussing the significance of the extreme lower tail is ihe 
fact that the time coordinate is given in Loq(t) and therefore 
4 small error of 0.01 order oi magnitude in Log (i) is 
translated in more than O.! order of magnitude error in useful! 
service life which in any ce is typically unacceptable. itum 
difficult to observe that thc exact determination of the extreme 
lower tail is based on how accurate one can determine Lhe Toad 
that could have been achieved If no failures паа пароепео соии 
the loading sequence. Therefore one can understand ла: ан 
cruisial importance to obtain exact data points with no irrelevant 
indications during the loading Seguence and Im 0700 
accomplisn this tas& two aspects have to De thoroughiv examined 


E 


г, 


How to avoid the irrelevant indications when evaluating 
the data points recorded during the loading sequence. 
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2 How to record as many failure data points as possicis 
«0 increase the precision of the final results. 


D STATISTICAL MODEL OF FAILURE 
The failure of a bundle composed of orittle fibers such as carbon and 
glass in a flexible matrix under tensile loads, is a complex statistical 


process involving the scattered failure of fibers at Impertection sites ana 


ол 
= 
~ 
(D 
Ce 
с) 
Cr 
Е 
ch 


Нэтс tne overloading of the adjacent 7ibers at tnese sites, a 
propagation of the neighboring fiber f3''ures to acritical size according to 
Rosen [ Ref. 6]. 

The rnathernatical model for the failure process was first presented 
by Rosen [ Ref. 6 | and is described in detail by Phoenix | Ref. 7 | However 
Гог purposes of completeness and in order to gain some insight in ihe 
features of the bundle failure a review of the simplest case is presentec 
herein. 

ее 21-50 2 ес 0 влее гпасегта! 15 considered to be 
Ee aO repeating pepresentative elements, the bundles, whicn 3l 
actually the focus of this work. Each oundle can conceptually be viewed as 2 


close planar (ог three dimensiona!) arrangement of n parallel fiber 
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filaments that form a tube as shown in Figure 6, which is loaded in 
tension. This tube structure is in turn conceptually partitioned into à ser!sg 


of m snort sections that are called microbundles. Furtermore the length 9 


of these short sections is called the meffective Jength after Rosen [Ref 


с 


6] and actually represents the length of the overload region of ir 
neighboring fiber filaments to a fiber failure site. In other words the 
composite consists of oundles each one of which has a length egual to mò 
and a volume of nmó. when a moderate tensile load per fiber p is applied co 
a bundle, fiber filaments fail randorly at sites where the strength is less 
than p and typically these random faliures are spatially quite tar appart. At 


the location of these failures sites the corresponding fiber filaments can 


not carry any load and therefore an overload of magnitude Kp takes place 


on each neighboring fiber filament, where K,>1 15 defined as in 


с) 


load concentation factor andr reoresents the number of consecutive 
failed fiber filaments immediately adjacent to a sound fiber Tieme 


counting on Doth sides. For this very simple case that is being сезер ы 


here Kr IS assumed to De: 
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INEFFECTIVE 5 2,» Злі? 
LENGTH | | 
РА т 
| FIBERS 
SULLY 
х STRESSED , 


| : tensile longitudinal stress 

foo : longitudinal fiber filament stress 

т 12: transverse in plane matrix snear stress 
longitudinal direction 


хааа 


Figure 6a Statistical Tensile Failure Mode! after ABosen,W.3. (Ref. 5] 





З : Ineffective Length of fiber filaments 
р : Load per fiber filament 
n Number of fiber filaments 
то: Number of fiber microbundLes 
О : FaiLure sites of fiber fiLaments 


Figure 50.  8und!e of fiber filaments in a matrix 
In tne form of a planar tape. rhe sundle 
1S COmoosed OF M Microounales in series 
eacn with n fiber filaments. Failure 
sites are locaiized within microoundles. 
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(ШЕ №. = 0/2 ОБ (Ба В. 


The length of this overload region which has already been defined as 8 15 


typically of the order of a few fiber filament diameters and both Кк. апа д 
depend on the fiber and the matrix material properties which can oe 
estimated by elasticity or numerical analysis. Under the increased load AnD 
additional failures will appear to the close neighboring fibers, which in 
turn will be subjected to a higher overload say Kgp (Kpp>K,p>i). This 


process continues until al! the failure occurences form a front of a critical 
size k*. lf this critical size is exceeded a catastrophic crack propagates 
throughout the bundle leading to the failure of the entire composite 
structure. 

In order to tailor this model to the thougnt experiment already 
described two cases of stress history of the composite are of interest. 


The first one is a constant stress history expressed as : 


CHES) o(t)=Q Гог 120 
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where Q is the constant positive stress level. This constant stress history 
is often referred to in the literature as stress-rupture of creep - 
rupture . Application of this type of stress history is schematically 
illustrated at Figure 4- 1a. 


The second one is a linearly increasing stress history expressed as: 


(11.6) со) є ВІ tor t20 


where R is a positive rate of stress increase. This type of stress history is 
assumed to account for the early failures of the fiber bundies that will take 
place during the loading sequence and it is often described as sfori - 
term strength. Application of this stress history can be viewed By looking 


at Figures 4-2a, 4-3a and 4-4a. The FCDF for the failure time of a single 
fiber filament of length 6 165 given bv Phoenix and Wu [ Ref. 4 ] after 


Coleman [Ref. 8] by the following relationship: 


(11.7) F(t;o) - -exp(-slo[ x(o(s)ds]) , for t20 


за 


where c(t] , 6 200 isthe stress history, 

кіх) = yx? |x 3 0 is the power-law breakdown rule 

o as well as 7 are positive constants, 

b{x] = axP , х 3 0 | із the Weibull shape function to 
compensate for the comrnonly observed Weibull — of the fibers “ШШ 
œ a5 uell as 8 are acain positive constants 

Experimental evidence [ Ref. 4}, suggests that 10 <p < 80 апа 2150 


that 1/4 < B < 4 for a wide range of fibers. It has also to be noted that & 


is proportional to the fiber-slement length, ог а = сб ior a Cons raft 


ag corresponding to the characteristic fiber length 8g. Considering д 10:06 
the characteristic length, since it is natural to the composite long single 
fiber filaments of length 1»» à willhave a taken as me where 
m  1/8.The theoretical justification of this fiber model 15 discussed ov 
Phoenix | Ref. 7 | in terms of the kinetic failure of idealized raoiecular 
Сгуо ии of the form found in stiff poivrneric filaments such as Kevlar riore 


specifically the key parameter p involves an approximation Го ine 
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potential function for chain scission and is shown to follow the 


~ 1 + 1. 
В nsnp: 


(18) юм Да 
where T is the absolute temperature, 
к 156 Boltzmann's сон... 


oe Pec oscars oi 2crivellon energy DUL IN ihe case 


а га 4 


of chain scission 1s Oniv about 3O 7$ of 
whereas the other parameters a , g and y vary with temperature in very 
complex ways. 

Essen theeesellangMNSEStress histories of interest described so 
far the behavior of single fiber filaments and fiber composites will be 
presented from both the sthrength and life standpoints together witn ine 
Dasic assumptions that are inherent to every case and finally some 
surintarizing results will also be presented to serve as a means of gaining 
Reece sca y Wotan 1OF TNE exper in cnlaiiy ankicipated results 3ccordina 


to Phoenix and Wu [| Ref. 4] and Phoenix: Ref. 5 | 


сл 
Cn 


|. Behavior of sinale fiber fitaments 
2. INO! CORT? DO 
The fiber filament strength is known to be closely 
approximated by the Weibull distribution which under the linearly 


increasing stress history [Eq.(11.6)] can be written in terms of the (FCDF) = 


ЕЖ as fallows: 


(119) ЕХ(р) = |-ехр(-р/р (ВУР(Р +?) "їе pO 


with a shape parameter g(o*1) and a scale pararneter equal to: 


(| За) рН) x m7 !/tB(p + 15 (В) 
where palR) = a РР“ (перу у АР pilot is ine scale 


- 


parameter for unit fiber length 5 апа Н 15 the loading rate. Finally the 


median fiber strength can be shown fo be equal to: 


56 


(11.9b) рж (В) - р(А)Пп(2)1 Ве“ 1)! 


Unfortunately the accuracy of Eq.(I11.9) in modelling the 


true FCDF can only be demonstrated for the lower tail of F*(p) that is for 
p««p, and this is because strengths associated with the middle and upper 


tail are not typically observed at laboratory gauge length | (of the order of 
several centimeters). However this does not constitute a serious problem 
since the results show that the knowledge of the middle and upper tails are 
not important in the characterization of composite strength. As it has 
already Deen explained what is really of practical importance is the 
accurate determination of the extreme lower tail. 


b. 5lress-hRupture 
A single fiber filament of length | = mô under the 


constant stress history described by Eq.(11.8) has a lifetime which can be 
described by the following formula after the appropriate reduction of 


БІН |17): 


(11.10) F(t) = 1-exp{-[t/t,(Q)}?) , for 120 


E 


with shape parameter g and scale parameter given by the formula: 


(11.102) 160) - m В (а) 


where %4(0) - о'0/(у2!'8) is the scale parameter for unit fiber 


filament length &. For the constant stress level Q the median time-to 


-failure can be evaluated as: 
(11. 10b) CRO « цох о(231178 


If the median strength Is Known at the stress rate R, and one wishes Lo 
calculate the median time to failure at the stress level Q, the equivalent 


median lifetime сап be evaluated by substituting Q in Eq (Il.10a) with 


ра (А) об Ед. (1.96) equal to: 


(ОГО) t*(R,p) = p% (R)[R(p+1)] 


зо 


2 Behavior of composite fiber bundles 


a. 2nort-term Sírength 


Опе сап now consider a single microbundle and define 
as  GQ(p) for p20 the (FCDF) of its strength or in other words the 


probability of failure of a stochastically selected microbundle under a 


nominal load per fiber filament equal (ор. Similarly one can also define 
Нар пера for p20 the FCDF for the strength of a bundle consisting of fiber 


filaments in a matrix. Since the necessary and sufficient condition for the 
survival of the bundle is the survival of each one of the rn rnicrobundies, 


the probability of survival for the whole bundle under the load p can be 
written as [1-G,(p)]™. Therefore the probability of failure of this bundle 


under the load np, amounts to the so called weakest link rule and can 


be presented as follows: 


(11.11) ШПЫЛ- І-ТОП | тог р>0 


Another basic assumption that has to be stated here is 


(he fact that each one of the surviving fiber filaments carries the same 
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amount of load a statement which amounts for the equa/-load sharing 
rule after Rosen [ Ref. 6 ]. The first ones to obtain accurate results for this 
kind of a model were Harlow & Phoenix [ Refs. 9 &1IO ] but only for small 
number of fiber filaments per bundle (i.e., n«14 ). Later the same authors 
[Refs. 11 & 12] developed a powerful recursion analysis for treatment of 
larger bundles. According to this analysis the FCDF of the bundle strength 
can be quite succesfully approximated by using the weakest link rule as 


follows: 


(11.12) Hg n(D) * 1-[1-W(QD)P | for p20 


where  W(p) is a very complicated function depending on the load 
concentration factor K. and the FCDF of the fiber filaments of length à 


under the tensile load p. The first one to develop a tractable approximation 
for W(p) was Smith { Ref. 13 ] who showed by means of asymptotic 
analysis for /arge n (ie mn = 109) that Eq. (11.12) can take the 


following form if one also replaces the constant stress history Q (described 
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in this case by the constant load p) by the linearly increasing stress 


history Rt: 


(1.13) жне (р) = 1-ехр(-[р/р (АК ВСР? 1)) , for p20 


with shape parameter  K*8(p*1) and scale parameter given by the 


formula: 


(1.132) рс (В) = (тп) КВО DI p (Ry dy 51 /(* D 


where K* represents the critical failure sequence size and mathernatically 


is the integer solving the equation : 


SD) I/r(k*1) « Bpg/In(mn) « 1/ r(X) 


where the function r(k) is defined by the following equation : 


(1130) г(К) 9 InC KK) - InC iL ПД ХК) апа 
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Гог Ко редици (11.134) 
AS before the median strength can be given by the equation: 


(1,14) р“ (R) = p. Gon c2)] / I BCo* D] 


D. Composite Stress-Hupture 
Similarly the FCDF for the lifetime of the composui 
material under the constant stress history described by Eq(ll.5) has а 


Weibull approximation of the form: 


(1.15) — *Hg y(t) = trexp(-[t/t (ay}k"F ‚ тог 120 
with shape parameter k*@ and scale parameter given by the equation: 


(11.16) 1 (0) = (та) КВЧ ot (0) 


where dyx has already been деѓіпеа бу Ед.(11.139) гог К = К апа ЕА(О) 15 
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defined in the same way as for the Eq.(Il.10a). The median lifetime of the 


model composite can therefore be given by the following expression: 


(11.17) ге (а) = (ап) КВ) 


Using the reasoning given in the analysis of the fiber 
filament behaviour if the median strength р“ „(В) is Known for the stress 
rate R the median lifetime under the stress level Q = р“ „(В) can be 


evaluated as: 
(1.18) t" C(R,p) - р (R/IR(G* 1] 


It has to be noted that the Weibull approximations (11.13) апа (11.15) 
are only valid when Bp is /arge, say greater than 6, as well as, that the 


best Weibull fit to Hm n(t) may produce a value of k* which is not an 


integer usually when the solution for k yields almost an equality on either 
side of Eq.(CIH. 13b). 
Summarizing the important characteristics of the model described 


one has to note that: 


|. Comparisons of the median lifetime and strength of a single fiber 
filament and a fiber composite are actually meaningless. 
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The variability in strength and lifetime for a composite is much 
less than that for a single fiber filament primarily due to the 
critical crack size. 


The size effect is a lot milder for the composite relative to that 


of the fiber for the strength and drastically lower for the lifetime 
again primarily due to the critical crack size. 
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НІ. THE CONTRIBUTION TO THE SOLUTION 


In accordance with the approach stated in the Introduction the 
required methodology appropiate to be used as a quide line in order to 
achieve the ultimate objective of obtaining strength-life data been 
presented in the of the problem approach section. 

For the equipment used by Carozzo [Ref. 2 : Fig. 3.4] shown in 
Figure 7 and from his evaluation of its performance (summarized in 
Figures B.! and B.2a of his thesis (IRef. 2]) shown in Figures 8a and 8b one 
can observe that although the desired target was the acquisition of data in ` 
hundredths of seconds (i.e 6000 data/min) the results suggested that a rate 
of only 3-4 data/min was achievable. In an effort to acquire more data 
during the loading phase, the testing traverse speed was lowered to the 
minimum. This slowest speed setting of the INSTRON testing equipment 
(at13 grams/sec) caused the drive motor of the INSTRON to stall erratically 
[Ref. 2 : pp. 29-30]. In addition, during the the creep rupture phase of the 
experiment the time dependence of creep was physically observed by the 


upward movement of the moving cross - head of the INSTRON. However 
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Figure 7  : Block Diagram of the initial Experimental Set-Up 
after Carozzo (Ref. 2] 
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Figure 8b : Load- Time curves for Test 22 
after Carozzo [Ref. 2]. 
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the data recorded to support this observation were considered completely 
unreliable in that it showed that the cross-head moved up and down which 
is equivalent to an entropy decrease ([Ref. 2 : p. 31]). These observations 
were what actually gave rise to the other goals of this work which are to 
evaluate the limitations of the testing equipment available and propose 
procedures that will serve in eliminating them (i.e explaining why the 
creep-rupture data were unreliable, avoiding the stall , by-pass the problem 
of inadequate accuracy of the extension measurements etc.), as well as to 
find out how more data could be acquired at least during the loading 
sequence for the reasons that have already been explained in the detailed 


analysis of the approach to the overall problem. 


A. TESTING EQUIPMENT LIMITATIONS 
The implementation of the INSTRON | model 1000 Universal Testing 
Equipment requires at least the description of it's operating principles 
prior to investigating it’s inherent limitations. 
ls Operating Principles 
The INSTRON 1000 is an 120 * 10V AC single phase 50-60Hz 


motor belt driven testing device capable of tensile and compressive 
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testing, extension measurements and the testing machine has been 
modified to perform creep monitoring if necessary. The load measurements 
are achieved by means of a combined Load Transducer (LT)/ Linear 
Variable Differential Transformer (LVDT) system, installed within the 
moving cross-head. The traverse crosshead measurement is achieved by a 
Shaft Encoder (SE) and the creep monitoring by means of a Digital 
Processor/Controller (DP/C) that can compensate for the cross-head 
adjustments relative to a predetermined load variation. The creep 
dispacement was specifically recorded by an additional LVDT that was 
mounted on the fixed cross-head of the INSTRON 1000. The operating 
principle of the LVDT is that of a transformer composed of a housing of 
the transformer coil and a core. The insertion of a ferrous rod-core alteres 
the electro-magnetic coupling thereby causing a change in the output 
voltage. Therefore by controlling the core insertion and measuring the 
output voltage one can evaluate a relationship that will allow for 
converting output voltage to displacement and vice versa. A detailed 
description of the major units, subsystems and enhancement features of 
the testing equipment and their primary functions is presented in Section | 


of Appendix A. 
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2 Determination of INSTRON 1000 Limitations by Simulation 


The strength of the fiber bundles to be tested is of the 50 105 
order of magnitude and the maximum utilization of INSTRON 1000 under 
the existing provisions set-up involves testing of nine (9) bundles 
simoultaneously (including the control bundle) Hence the load level of 
interest is of the 450 Ibs order of magnitude, an estimate that dictates 
Гог the use of only the 1000 Ibs LT during the simulation. Given also that 
the graphite fiber bundles ought to be tested in tension, because the area 
of interest is ultimately the broad class of composite applications where 
the structure is Sub jected to sustained tensile loading, as well as that the 
graphite fiber bundles are expensive to use them for purposes other than 
the actual experiment a simulation process is obviously necessary. During 
this simulation experiment, the above mentioned load level (450 Ibs) will 
have to be easily achieved and/or overcome a good number of times so that 
several conclusions could be drawn as far as the limitations of the 
equipment are concerned. 

The objectives, the theoretical model, the design of a simple 


mechanical set-up that had to be used with the INSTRON 1000 during the 
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simulation and the detailed procedure of this aside task are presented in 
Section Il of Appendix A. 
3. Stability of the LVDT 

The background for checking the stability of the LVDT lies on 
the fact, that although the cross-head of the INSTRON must have moved 
during the creep-rupture phase, no reasonable output had been recorded 
accordingly and the output (displacement) was measured by means of an 
LVDT which was externally mounted on the INSTRON 1000. 

The background, the objectives, the set up used and the 
detailed procedure of this separate aside task are presented in Section |!!! 


of Appendix A. 


B. DATA ACQUISITION SYSTEM LIMITATIONS 

One of the major problems associated with the graphite veni life 
testing as presented by Carozzo [Ref. 2: Fig. B1 & B2a] (Figures 8a and 8b) is 
the slow data aquisition rate that was achieved during the loading sequence 
(0.04 readings/sec) was substantially less than the manufacturer's 
specification of about 2 readings/sec [Ref. 15 : p. 32]. For completeness and 


future reference purposes the experimental set-up used by Carozzo during 
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the initial run of the experiment is prcsented in Figure 9. A survey of the 
available literature on this subject revealed that the reading rates can be 
influenced by many factors. 
|. Signal environment factors 
The most significant factors in this area are the line related 
factors (the smaller the loop the better , 60 Hz line frequency gives about 
172 faster reading rates), [Ref. 15 : pp. 30,34] & [Ref. 16 : p. 42], the 
broadband noise (the larger the loop, the more the noise contrbutions from 
the various accesories in the loop), [Ref. 15 :p. 30] & [Ref. 16 : p. 42] and the 
thermal gradients (negligible in this case due to the relatively constant 
temperature of the lab environment), [Ref. 15 : p. 30]. 
2. Desired accuracy 
At 3.5 digits of resolution ( N3) the ability of the voltmeter to 
acurately measure DC voltages in the presence of AC voltages at power line 
frequencies is expressed by the Normal Mode Rejection (NMR), [Ref. 15 : pp. 
30, 43], which in this case is O db NMR. The smaller the number of digits 
of resolution the smaller the "Integration Time" which obviously determines 


the reading rate and hence for (N3) the reading rate is faster. Similarly at 
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Figure 9 : Block diagram of the experimental set - up 
used during the initial run of the experiment 
Involving one INSTRON 1000 machine. 
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N4 (59 db NMR) ап sntermediate reading rate is achieved and finally at NS 


(80 db NMR) the reding rate is s/ow. 


(1) 


(2) 


(1) 


(2) 


(3) 


Ч 


Type of measurement 


a. Resistance measurements [Ref.15 :pp. 30-32] 


At 30 Mohm ( 200 msec settling time ) the recording rate is 1.5 
readings/sec. 


At 3 Mohm ( 20 msec settling time ) the recording rate is 2.1 
readings/sec . 


D. Voltage measurements. 


The DC readings are S to 50 times faster than the AC ones [Ref. 15 
ED. 32 


Positive voltage measurements increase the reading rate. Possible 
problems can be raised by the LVDTs whose outputs can be 
negative dependent on their null point with respect to the 
measurements under evaluation. 


AUTOZERO is a function that allows the operator to enable or 
disable the internal zero correction of the voltmeter of the 
HP-35421A. The default condition is ON ( used by the set-up of 
Figure 9 ) but turning AUTOZERO OFF could substantially 
increasthe reading rate but with a tradeoff of long term stability 
and accuracy [Ref. 15 : p. 41]. However no problems of this sort are 
of interest during the relatively short period of the loading 
sequence which is of the order of a few minutes (see Figures 8a 
and 8b). 
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4. Inherent to the HP-3421A features 

Calibration has to be performed after the installation of the 
one 44462A / 10 channel Multiplexer (option 020) but it is of interest only 
when absolute values are necessary; however this is not this case. 

AUTORANGING is a function meaning that signals with a wide 
dynamic range such as switching from channel to channel to take readings 
are made more quickly by utilizing it [Ref. 15 : p. 44]. But this function is 
made primarily for taking readings in the N4 mode and not in the N5 mode 
which has to be used due to the HP-41CX involvement(see Figure 9). 

When operating with the HEWLETT PACKARD-INTERFACE LOOP 
(HP-IL) most 3421A commands will " Hold - Up " the computer until all 
readings have been taken. These problems are known as information 
transfer ones. Recording time by the HP-41CX when the order of magnitude 
of the 3421A capacity is about 2 readings / sec and the target is to take say 
1 00 readings using one 444624 10 channel multiplexer assembly 50 seconds 
are required for the 100 readings. The computer will not be able to perform 
any operation during the 50 seconds while the 3421A is taking the 


measurements. 


76 





The 3421A commands that do not hold-up the computer are the 
so called advanced commands of digital monitoring (MN, MH and ML) the 
totalize function (TOT) and the digital trigger (DT) commands. However if 
any of the digital monitor modes (MN ,MH ,ML or DT) are in effect in the 
3421A any communication through the HP-IL will cause the mode to be 
aborted. In other words if the 5421A is waiting for a digital trigger and the 
computer or another instrument in the loop sends any commands or data 
through the interface, the 3421A will abort the digital monitor and will no 
longer respond. 

On the other hand, according to Carozzo [Ref. 2 : p.25] the basic 
data logger software program provided by the the HP-44468A Data 
Acquisition Control Pac was modified to permit the recording of time in 
hundredths of seconds and edited to reduce the number of data registers 
uscd. However recording time in hundredths of seconds is 50 times faster 
the maximum capacity of the HP-3421A DACU (2 readings/second). In 
addition to this the system timing is greatly dependent on the nature of the 
set-up itself and especially on whether the printer or tape cassette (or 
both) are used the number and type of measurement sequences, User 


Functions as well as the number of channels in each sequence [Ref. 16 :p. 42]. 
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Several contacts with system engineers of Hewlett-Packard revealed that 
the reading rates of the HP-3421A manual for the case of the experiment 
under evaluation, (21F1-2.19/13.18/22.22) [Ref. 15 : p. 32] refer only to the 
unit itself and if this unit is interfaced with the HP-4l calculator they 
reduce to 1.0/1.0/20 whereas no data were available in case the printer 
and the cassette drive are within the loop. Their suggestion with respect to 
the very slow data acquisition rate was to definitely take out of the loop 
the printer and the cassette drive and also try to use AUTOZERO OFF and 
AUTORANGING, or switch to another system as the HP-85 desk top computer 
or the neuer HP-71 which provide the opportunity of using the advanced 
commands that allow for better programming flexibility (ie. use of 
AUTOZERO and AUTORANGING functions and as less digits of resolution as 
possible according to the requirements of the experiment) and also by-pass 
the problem of storing and printing data by using their own modules 
without contributing to a larger interface loop. Under the above 
circumstances (interfacing the HP-3421A DACU with the HP-85 and using 
the advanced commands) according to the HEWLETT-PACKARD engineers the 
recording rates of the HP-3421A DACU as presented in the Operating, 


Programming and Configuration Manual [Ref. 15 : p. 32] reduce to 3/9/11 


78 


respectively depending on the number of digits of resolution required which 
Is substantially higher to what already has been achieved. 

The points of contact with HEWLETT-PACKARD as well as the 
official correspondence made with the customers division, by which the 
company was asked to confirm by letter the information released by phone 


and suggest possible solutions of improvement is presented in Appendix B. 


t OVERALL EXPERIMENTAL SET-UP 

To come up with the overall experimental set-up several caclulations 
had to be made concerning the optimum, yet safe number of pieces of 
equipment that ought to be used. These calculations are presented in 
Appendix C and dictated the utilization of one (1) power supply НР-62168 
for every set of nine (9) Load Cells ( per INSTRON 1000 ) as well as an extra 
power supply HP-6216B for the LVDTs that would be externally mounted to 
the INSTRON machines for the recording of the creep-rupture displacement. 

The basic steps that were followed for the overall experimental 
set-up which is presented in a diagram form by Figures 10a and 10b, are 
also given in Appendix C under the title " Implementation of the Overall 


Experimental Set-Up ~”. 
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Figure 104: Block diagram of the overall experimental 
set-up involving four INSTRON 1000 machines 
oroposed for the Loading Sequence phase. 
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Figure 106 : Block diagram of the overall exper'mental 
set-up involving four INSTRON 1000 machines 
proposed for the Creep-Rupture phase. 
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А. INSTRON LIMITATIONS 

According to the INSTRON manual [ Ref. 14 ] the Liquid Crystal Display 
(LCD) indicator for the extension readout of the moving cross-head of the 
INSTRON between the no-load case and the under load situation, gives an 
agreement within about 10 % and hence only the latter, which is directly 
related to the ultimate experiment, was thoroughly examined. 

Several attempts were made to obtain displacements for the 
minimum knob speed setting (i.e 0.5 in/min-MIN) but it turned out that no 
movement of the cross-head had taken place as it was observed by both the 
LCD and the DIAL indicators; the latter was used to monitor the physical 
displacement as described in the Appendix A. Furthermore it - observed 
that after a short period of time the LCD displacement readout started to 
decrease with time. 

Using the bracketing technique it was found that no matter what the 
load was (within the capability of the equipment) no stall occured for knob 


speed settings larger than 0.5 in/min-30 % and the overtravel-limiters 
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system had finally to interrupt the loading sequence to protect the 
equipment from overload. 

A good number of attempts were also performed at the range of 
speeds (0.5 in/min-MIN, 05 in/min-12 2) This showed completely 
irrelevant indications (Table 1) as far as the observed and the expected 
spceds correlation was concerned and suggest that in actuallity the 
cross-head had moved so little and in such a random way that no 
Significance could be given to these small Variable Speed Indicator (VS!) 
settings. 

several knob speed settings were tried to evaluate the quality of the 
response and some representative sample data are given in Table I! which 
actually turned out to be in very good agreement with the anticipated 


results. One can also observe from Figures 11a, 116 and llc 


corresponding to the above mentioned table that the functions d pja = КО 


d, cp * f(t) and P » f(t) are single valued, monotonically increasing and as 


the time increases there is a tendency for the displacement and the load 
intervals to decrease. Having thus verified the anticipated results and 


observing from the quality of the response, that the helical spring employed 


83 


TABLE | 


INSTRON 1000 RESPONSE FOR LOW SETTINGS 
OF THE VARIABLE SPEED INDICATOR 


Column 1 Сота 2 Coiumn 3 Column 4 Сойлтт 5 Column & Coiumn 7 
2:2 0.025 224 12.0 61 1 125 
28 0.040 254 1157 61 Ї 20.0 
5.6 0 060 194 ШЕЛ 61 1 30.0 
20 0.075 254 11.7 61 1 325 
3.0 0.090 591 30.0 63 0 45.0 
3.6 0.100 50 3.0 61 2 50.0 
1.4 0.110 168 4.0 64 5 55.0 

Ї 25 0.110 167 315 63 2 55.0 
9.7 0.110 118 19.0 62 2 55.0 
12.6 0.110 76 16.0 62 1 55.0 
8.9 0.110 945 140.0 85 15 50.0 


CBSRVD SPO VREL SPO IND TOTAL TIPE OAL OSPLCT STL LO (P) LCD OSPL EXPCTD SPD 
x(.001in/min) DIAL SETTING x ( sec ) х(.001 а) x( ibs ) x(.01 i» ) x(.OOlin/min) 
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ASCE || 


TYPICAL DATA BASE SHOWING THE QUALITY OF THE RESPONSE 
OF THE INSTRON 1000 TESTING EQUIPMENT 


Column | Column 2 Column J Column 4 
ЛЕ 0 50 0 
29 4 61 0 
50 9 61 1 
.75 15 62 1 

1.00 22 55 2 
125 27 63 3 
1.50 32 54 3 
1.75 57 65 4 
2.00 44 66 4 
272 50 67 З 
2.50 56 68 5 
2:75 62 69 6 
2.00 58 71 7 
3225 79 71 7 
3.50 80 73 8 
3.75 64 73 9 
4.00 86 74 9 
423 92 75 [0 
450 98 77 10 
4.75 104 78 11 
5.00 119 79 12 
dui 116 80 12 
5.50 122 81 15 
3.15 129 53 14 
6.00 134 84 14 
625 159 65 13 
6.50 145 87 14 
6.75 150 88 14 
7.09 155 89 15 
723 162 91 15 
750 168 92 16 
775 174 93 16 
8.00 180 Ж 17 
825 182 24 18 
8.50 187 95 18 
8.75 193 97 19 
9.00 195 397 16 


TIPE OAL OSPLNT ТУ LOAD (СО OSPLHT 
x (minutes) x(0.001 incn) x ( lbs ) x (0.01 inch ) 


COOES : *: STAL LOAD 
DATA IDFN : RUM * 18 
YSI STN6 :  Z2x.13 in/min 
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lla: 


з ~ 0.5726 + 22.254x R= 1.00 


RUN *18 ( VSI : 0.5x0.13 in/min ) 





Time (min)-Column 1 


Typical response for the Loading-Sequence. 
Variation of physical displacement (DIAL) 
of the INSTRON 1000 moving cross - head 
with time. 
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у =0.1598 + 2.143х В=0.99 


О Column 4 


LCD DspImt (0.01 in) 


RUN © I8 ( VSI : 0.5x0.13 in/min ) 





Time (min)-Column 1 


Figure llb: Typical response for the Loading-Sequence. 
Variation of displayed displacement (LCD) 
of the INSTRON 1000 moving cross - head 
with time. 
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Tensile Load (lbs) 


RUN * I8 (VSI : 0.5X0.13 in/min ) 





Time (min)-Columa 1 


Figure llc: Typical response for the Loading-Sequence. 
Variation of tensile load applied to the 
INSTRON 1000 moving cross - head 
with time. 
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managed to represent quite well the linear model assumed, the data base of 
Table Il! was created by repetition of the loading sequence at various 
speed settings. A closer look to this data base yielded the following results. 

There is always an one to one correspondence between the following 


pairs of variables of interest: 


|. The LCD and the DIAL displacements. Their relationship appears 
to be a linear under a slope of one (1) as shown in Figure 12a. The 
greater error (noise) band which is observed to the lower side of 
the simple curve fit is attributed to the sequence under which the 
data were manually recorded. Had the data been recorded in the 
opposite sequence, a greater error band to the upper side of the 
simple curve fit would have heen observed. 


2. The expected speed and the observed speed as shown in Figure 
120. 

57 The stall load and the Variable Speed Indicator (VSI) as shown in 
Figure 12с. 


The plot of stall load versus the VSI as presented in Figure 12C 
enables the descission maker to predict the appropriate knob combinations 
to reach a certain load level without stalling. However this plot is 
characterized by a lot of scatter in the data and that is why it is also 


presented in a form to compensate for the mean upper and lower values 


89 


TABLE III 


DATA BASE OF THE STALL CONDITIONS DURING THE SIMULATION 
FOR VARIOUS SETTINGS OF THE VARIABLE SPEED INDICATOR 


Column | Calumn 2 Column J Column 4 Colum 5 Column 6 Column 7 
9.13 55 20 15 13 аку 34 
0.13 65 146 197 42 799 97 
9.13 55 193 638 55 1961 201 
9.13 65 22.9 195 18 552 97 
0.13 55 19.4 21 0 55 52 
9.13 65 223 71 6 17 70 
9713 55 12.9 55 3 245 69 
0.13 63 13.0 55 5 254 69 
0.15 55 12.9 56 3 250 69 
0.15 65 14.1 53 5 234 69 
9.13 65 22.4 205 19 550 37 
9.14 70 258 749 TS 1223 214 
9.14 7 29.5 411 41 834 148 
9.14 7 252 194 18 età: 95 
0.14 79 720.5 165 13 323 93 
0.14 7 513 SQ 5 141 68 
9.14 70 40.5 57 167 89 
9.14 70 22 90 7 243 76 
9.14 70 199 147 15 445 as 
9.14 79 238 37 4 144 99 
0.14 70 213 110 8 303 90 
0.14 ' 70 237 69 5 122 99 
9.14 70 18.1 155 13 S47 92 
9.15 d 232 973 79 2404 279 
9.1$ 75 34.0 425 43 750 152 
0715 T 15.4 1125 112 1925 321 
0.15 25 217 366 87 1638 250 
0.15 io 20.7 1125 197 2519 T 
9.13 75 33.0 1200 125 2185 31821 
0.13 75 345 911 35 1555 259 
0.13 79 251 1529 143 2815 423 
9.15 ТЗ 332 173 13 235 90 
9.13 7S 53.9 742 70 1505 215 
9.13 73 2153 182 18 400 3-4 
0.15 75 58 112 11 251 81 
0.15 79 245 168 14 414 92 
0.15 25 242 394 36 977 142 
9.16 20 53 1559 153 60 456 
0.15 50 33.9 1557 149 2449 424 


УЫ STNG ЕХРСТО SFO OBSRVUSPO CAL OSPUTT (0 САТ ТОТА ITE STALL LOAD 
x (0.5 питта) x(.OOlin/rmm)  {d.00 liw ma) z (001: ) х ( 91 т; x (sec ) х (1) 
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TABLE 11 


DATA BASE OF THE STALL CONDITIONS DURING THE SIMULATION 
FOR VARIOUS SETTINGS OF THE VARIABLE SPEED INDICATOR 


Column | Column 2 Сацята 5 Catumn 4 Coiumn S Column 6 Сойятт 7 
0.16 50 397 960 90 1451 270 
0.16 20 32. 50 5 53 62 
9.16 80 42.0 963 101 1375 274 
9.16 89 4.1 165 16 225 30 
9.16 30 23-31 287 +4 822 166 
0.16 9% ща SS! 59 993 185 
0.16 до 24.4 770 82 1344 
9.16 до 5.9 569 81 1537 1-4 
0.16 80 ша аф 41 2193 472 
0.16 50 239 716 64 1258 2! 
9.17 65 497 1467 150 1770 
0.17 65 423 1796 159 2316 437 
9.17 85 43.7 226 93 1165 245 
9.17 83 465 414 41 5-6 143 
9.17 85 42.7 1542 120 1255 156 
9.17 65 «2 410 81 310 193 
0.17 T 415 717 74 990 <0 
0.17 65 43 3 779 73 1054 в 
0.17 23 40.4 995 94 1477 234 
0.17 55 408 1020 33 1543 239 
0.17 85 372 1214 92 1954 о 
0.18 90 47.5 1681 176 2126 266 
9.18 90 49.9 1695 160 2037 467 
0.18 99 = 5 1155 ГІЗ 1287 247 
0.18 90 $19 1441 138 1566 400 
0.18 90 20.7 1220 73 1960 379 
0.18 90 539 989 98 1104 284 
9.18 90 49 3 1272 117 14357 3140 
0.18 90 44 0 1352 9? 1513 759 
0.18 90 49 8 767 72 775 224 
9.18 90 29.5 906 87 927 2 
9.15 90 49 3 1346 151 1632 za 
9.18 90 S23 1210 115 1557 340 
0.19 95 69.3 515 35 507 172 
0.19 95 779 2021 304 1557 552 
9.19 95 25.5 2195 187 2277 592 
0.19 95 482 1651 174 4057 462 
9.19 25 23.5 624 81 627 242 
9.19 95 572 915 за 957 282 
9.19 95 33 4445 6! 873 200 
0.19 95 46 6 1775 257 2197 489 


VSI STNS КОСУ 5РО CBSRVD SPO ДА. ОА МТ LOD OSAT TOTAL ПРЕ STALL LOAD 
x (0.5 in/rmn) x(.OOtin/rmn) — х.ОО Ппитва) x (001: ) «( Otm) 1 (sec ) х (№) 


М SS eee et ee Se ee ЫЕЕЕ 
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TABLE iii 


DATA BASE OF THE STALL CONDITIONS DURING THE SIMULATION 
FOR VARIOUS SETTINGS OF THE VARIABLE SPEED INDICATOR 


Column 1 Column 2 Column J Column 4 Column 5 Соната б Column 7 
0.19 95 43.1 1577 145 1537 482 
9.19 95 433 1566 ics 2474 513 
0.19 95 575 390 37 307 141 
0.19 95 70.1 1274 107 1091 250 
0.19 95 47.8 65 10 32 59 
020 100 973 543 % 747 248 
020 100 69.5 1250 129 1077 344 
020 109 655 1453 124 1317 571 
0.20 100 323$ 2005 241 3224 608 
920 100 70.1 157 15 143 91 
020 100 50.2 1152 113 1148 323 
020 109 478 1306 50 1647 259 
029 109 56.9 1050 114 1197 295 
0.29 100 733 444) 25 579 149 
020 100 511 2145 183 2519 553 
0.29 109 47 9 16-43 199 4515 324 
9.20 100 50.1 1600 130 1917 «357 
029 100 S 770 72 777 249 
029 109 49 5 1556 113 1287 423 
021 105 TX тео 213 1227 912 
021 105 685 1048 115 918 306 
0.21 ics 59.4 1642 153 1659 451 
021 105 633 1359 113 1270 769 
921 105 49 1710 158 2097 474 
0.21 195 58.4 1701 252 1747 474 
221 105 412 1965 129 1363 545 
2:22 110 75.9 2 5 57 70 
0.2 110 56.1 1927 169 1658 249 
9.22 110 30.6 321 92 511 240 
0:22 110 1015 1813 | 44 1972 400 
222 110 773 1149 111 536 322 
0.23 115 742 2270 ara 1917 TS 
0.43 115 99 9 1865 181 1120 514 
023 115 69.7 2116 193 1822 596 
0.23 115 36 В 495 50 30 166 
0.25 tis 515 2227 50 423 168 
0.24 120 915 2240 216 1468 783 
924 120 90.4 545 58 S26 453 
924 120 109.7 425 40 233 125 
0.24 129 71.4 1410 140 1165 331 


VS} STNG ОФСТО ЗРО  OBSRVOSPO МА, О5РЈТП (ОО ОРТ TOTAL TTE STALL LOAD 
x (0.5 m/man) x(.OOlie/mem)  x( OO lin/mmn) x (001! ) (С Olin) x (sec ) x (Ms) 


м шї] тағ По ~ НЕ Шо EE алаңын NNN о НННП 
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550 


у = -1.1254 + 0.0979х R=0.95 


500 


250 


200 


150 


100 


LCD ОЗРЕМТ ( O.01 іп) 





Q 300 1000 1300 2000 2900 3000 3500 


DIAL DSPLHT (0.001 in)-Column 4 


Figure 12a : Quality of the Response of the Simulated Experiment. 
Relation between physical (DIAL) and displayed (LCD) 
displacement of INSTRON 1000 moving cross-head at 
stall conditions. 
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OBSRVD SPD (0.001 in/min) 


Figure 





125: 


= - 67.4375 + 1 2851x R* 0l 


Column 3 


70 80 90 100 110 120 130 


EXPCTD SPD (0.001 in /min)-Column 2 


Quality of the Response of the Simulated Experiment. 
Relation between observed speed (OBSRVO SPD) and 
expected speed (LXPCTD SPD) of INSTRON 1000 
moving cross-head at stall conditions. 
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800 


у = 20.222 * 107(5.9459x) В з 0.60 


600 <> Column 7 


400 


STALL LOAD (1^5 x) 


200 





UNE 0.14 0.16 0.18 020 0222 0.24 


VSI KNB STNG Ө O03 in/min-Columna 1 


Figure I2c : Quality of the Response of the Simulated Experiment. 
Relation between achieved stall load and INSTRON 
|000 variable speed indicator knob setting (VSI KNB 
STNG). 


са 


values achieved at every VSI knob setting as Figure 12c-1 after the 
appropriate data manipulation presented in Table III-A. It is therefore 
evident that in order to avoid stall one has to take into account Ше VSI 
knob setting corresponding to the mean which also gives the higher 
confidence level in that no stall will occur. For the actual experiment 
where nine (9) fiber bundles will have to be tested simoultaneously totaling 
450 lbs order of magnitude load, Figure 12c-1 suggests that the VSI knob 
setting to be used has to be that of 0.5 in/min-21 9. 

The graphs of Figures 12a, b, c, c- | are characterized by a great 
amount of scatter which appears to be of increasing magnitude as the VSI 
Setting is increased. However this is not true and this indication is only 
attributed to the limitations of the testing equipment itself. Had the 
Capacity of the INSTRON been larger the scatter would decrease almost to 
zero. This statement ows its credit to the large amounts of attempts that 
have been made for VSI knob settings larger than 0.24 in which the switch 
limiter-system had always taken over to interrupt the loading sequence for 
safety purposes in the neighborhood of 950-975 Ibs. Being conservative this 
means that a stall load of at least 1000 Ibs for the VSI of 0.25 and 1200 


lbs for the VSI of 0.26 would have been achieved if capacity were not a 
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TABLE III-A 


MIN, MEAN & MAX VALUES OF PARAMETERS OF INTEREST 
FOR VARIOUS VSI SETTINGS AT STALL CONDITIONS 


Column 1 Column 2 Column 3 Column 4 Column 5 
' 6.15 65 2:0 16.0 223 
0.14 70 18.1 19.6 305 
0.15 79 257 30.7 55.4 
0.16 80 4.1 315 42.0 
0.17 85 372 435 49.7 
0.18 90 40.7 2022 58.6 
0.19 95 43.1 60.6 77.9 
0.20 100 325 58.0 78.5 
0.21 105 41.2 58.9 725 
0:22 110 75.0 502 101.5 
0.23 118 69.7 82.8 99.9 
0.24 120 71.4 90.8 109.7 


VSI STNG EXPCTD SPD MIN OBS SPO MEAN CBS SD MAX CBS PSD 
(0.5 in./min) x(.OOlin/min) x(.OOli/min) x(.OOtin/min)  x(.OO lin/min) 
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MIN, MEAN & MAX VALUES OF PARAMETERS OF INTEREST 


TABLE 


111-А 


FOR VARIOUS VSI SETTINGS AT STALL CONDITIONS 





Column 6 Column 7 Column 8 Column 9 Column 10 Column 11 
150 141.9 638.0 0.0 14.5 65.0 

50.0 107.8 411.0 4.0 АТ 72.0 

20 714.1 1929.0 110 67.4 145.0 

50.0 759.4 1587.0 3.0 72.B 1530 

414.0 10221 1796.0 41.0 97.7 159.0 
7579 12552 1695.0 720 116.6 176.0 
65.0 1209.3 21050 10.0 130.8 304.0 

167.0 1258.1 2145.0 15.0 176.6 941.0 
1048.0 1658.9 22070 115.0 161.9 25240 
79.0 112728 1827.0 5.0 106.0 188.0 
495.0 1486.6 25700 50 142.8 2250 
426.0 1259.8 2240.0 40.0 12180 216.0 

MIN МА MEAN DIAL MAX DIAL те 100 MEAN LCD МАХ LCO 
x (.001 in) x (.001 in) x (.001 іп) х (01110) х (.01 іп) x -(01: 2 
CONTINUED PAGE 2 CF $ | PAGES 
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TABLE 


LISA 


MIN, MEAN & MAX VALUES OF PARAMETERS OF INTEREST 


Column 13 


503.4 

465 5 
14395 
1555.6 
14275 
1929.7 
1558.6 


Column 14 


1961.0 
1555.0 
2515 0 
2560.0 
256,0 
2126.0 
2474.0 
3224.0 
2865.0 
1658.0 
1917.0 


Column 15 


62.0 
69.0 
81.0 
62.0 
145.0 
224.0 
69.0 


FOR VARIOUS VSI SETTINGS AT STALL CONDITIONS 


Column 16 Colurnn 17 


89.5 
105.1 
2193 
240.7 
295.2 
552.5 
505.4 
541.4 
461.3 


201.0 


MIN TIME MEAN TIPE MAX TIME MIN STL LD MEAN STL LD MAX STL LD 
x ( sec ) x ( 15 ) x ( Ibs ) х ( 105) 
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y = 77.9514 * 10%4.2009х) Я = 0.86 
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Figure 12с-1 : Response Quality of the Simulated Experiment. 


Relations between MIN, MEAN and MAX achieved stall 
load and INSTRON 1000 variable speed indicator 
knob setting (VSI KNB STNG). 
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problem. Including these postulated values in Table III-A one can come up 
with Figure 12c-2 which clearly shows the decrease of the error band 
mentioned above. 

The accuracy of the LCD is limited to 0.0! in and therefore 
displacements of lower order of magnitude cannot be correctly recorded as 
for example during the graphite fibers creep. Roughly the deformation 
during creep is about 20% of that achieved at the end of the elastic 
behaviour. Knowing that the elastic deformation is of the order of magnitude 
of 2 %, for the gauge length of 10 inches to be used, the ultimate creep 
displacement would be of the order of magnitude of 0.04 in and therefore it 
is unlikely that the INSTRON LCD could record any accurate data during the 
creep experiment. 

Another point worthy of comment is that from the simulation of the 
actual experiment via the helical spring, it was observed that whenever a 
stall load level was reached the LCD display appeared to decrease in time, 
the order of magnitude of which was estimated to be about 0.01 іп / 20 sec. 
Hence the LCD display cannot be considered a reliable means of recording 
the fiber's creep data which is of primary importance for life testing. The 


reason for this decrease in the LCD displacement reading lies on the fact 
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Response Quality of the Simulated Experiment. 
Convergence of stall load relations with VSI KNB 
STNG under theoretically extrapolated conditions. 
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that this reading ís achieved through the Shaft Encoder (SE), a device that 
generates an output in proportion to the revolutions per minute of the 
Intermediate Drive Shaft (IDS) by which it is driven. Consequently in case 
of a stall, due to the internal friction, the IDS is not rotating, no output is 
generated ру the SE and hence the decrease in the LCD displacement 
reading. 

Finally summarizing the validity of these results one can rely on the 
INSTRON 1000 only for the loading sequence until the stall load is achieved 
by merely selecting VSI knob settings as dictated by the relatioship of 
Figure 12c-1, but cannot rely on the LCD for the creep-rupture phase and 
QT other means of recording and/or displaying the displacement 
during the creep experiment have to be implemented. Such means is the 
SCHAEVITZ LVDT (for recording) mounted on the loading frame of the 
INSTRON in conjunction with a high accuracy multimeter as the available 


FLUKE 8440A to measure and display it's output. 


D. LVDT LIMITATIONS 
The bracketing technique outlined in Section II! of the Appendix A 


under the title “Procedure” showed that the LVDT was totally insensitive 


103 


for an input of less than 11.0635 V. This result means that for al! practical 
purposes the зай is too low for being able to be recorded and/or 
displayed for an input less than the 12 V order of magnitude. On this basis 
the LVDT response was then checked at three points namely 24 V, 18 V and 
12 V respectively, for the. displacements of 0.02 in, 0.04 in and 0.06 inches 
due to the fact that during == the order of magnitude of the creep 
displacement for the graphite fiber bundles is known to be of the 0.04 order 
of magnitude ( 2 % of 20 & of 10 in ). 

For a positive direction of displacement ( i.e 0.02, 0.04, 0.06 in ) and 


after averaging the output of about 10 readings per case, one can come up 
with Table IV-A and observe that , the closer the input voltage (Vi,) to 


the calibration value of the 24 V, tne better the response and vise versa. In 
fact the normalized response of the LVDT at the end points of the voltage 
region (18V-24V) showed almost exactly the same slope, Ж it was 
anticipated and a high goodness of fit as shown in Figures 135a and 13b. 
For a negative direction of displacement ( i.e 0.06, 0.04, 0.02 ) one 
can similarly come up with Table IV-B which leads to exactly a 


similar observation as one can see by looking at Figures 13c and 13d. 
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TABLE IV-A 


LVDT OUTPUT FOR VARIOUS INPUTS AND POSITIVE 
DIRECTION OF DISPLACEMENT 


Column 1 Column 2 Column 3 Colurnn 4 Column 5 Column 6 Column 7 
0.02 -1.76213 -0.14684 -3.52404 -0.19578 -5.57715 -0.23238 
0.04 -1.71797 -0.14516 -3.54925 -0.19218 -5.46213 -0.22759 
0.06 -1.31844 -0.10987 -3.29234 -0.18291 -5.27860 -9.21786 


OSPLCI"NT(in) AVRG Vout(V) Vout / Vin AVRGVout(V) Vout / Vin AVRGVouti V) V out / Vin 
Pstv Orctn Vin z 12V Vin » 12V Vin = 18V Vin » 18V Vin » 24V Vin s 24V 


 —Ó——— — ны ————— € ——  ——— |  ——M— «s A Ee eee 


TABLE IV-B 


LVDT OUTPUT FOR VARIOUS INPUTS AND NEGATIVE 
DIRECTION OF DISPLACEMENT 


Column 1 Column 2 Column 3 Column 4 Column 5 Column 6 Column 7 

0.06 -1.63297 -0.13508 -3.30045 -0.18336 -5 24423 -0.218S1 

0.04 -1.71655 -0.14505 -3.45410 0.19189 -5 46877 -0.22787 

0.02 -1.76616 0.14718 -3.53410 -0.19634 -5.59618 -9.23317 
OSPLCIINT(in) | AVRG Vout(V) Vout / Vin AVRGVout(V) Vout / Vin AVR6 Vout(V) V out / Vin 
Ngtv Orctn Vin a 12V Vin 3 12V Vina |8М Vin a 18V Vin » 24V Vin z 24V 


SSS LS я 0 ee 
= 
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Average Veut (V) 
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у= - 2.0432 + 11.0923x В= 0.91 






(Vin = 18V) 


y = - 3.6202 + 5.7925х В = 0.96 
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* Column 4 
-$ и Соџтпбо 
(Vin » 24V) 
y = -5.7712 +8.7137x R=0.98 
-6 
0.01 0.02 0.03 0.04 0.05 0.06 


Pstv Displacement (in}-Columa 1 


Figure 13а: Response of the УОТ for various inputs 
and directions of displacement 
Variation of Average Vout with Positive 
displacement. 
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GO Column 3 
¢ Column 5 (Vin = 12V) 


и Column 7 


70.15 у = -0.1703 + 0.9243х 8=0.91 


(Vin * 18V) 


-0.20 у= -0.2032 + 0.3217x R=0.97 


(Vin = 24V) 





у= - 0.2405 +0.363x R= 0.98 


0.01 0.02 0.03 0.04 0.05 0.06 
Реб Displacement (in}-Column 1 
Figure 13b: Response of the LVOT for various inputs 
and directions of displacement. 


Vartation of ( Vout / Vin) with Positive 
displacement. 
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-2 y = - 1.8384 + 3.3298х В = 0.99 


Average Vout (V) 





-3 
(Vin » 18V) 
у * -3.6632 + 5.8412х В*0.98 

-4 

2 Column 2 

¢ Column 4 

8 Column 6 
"à (Vin - 24V) 

у = - 5.7883 + 8.7988x В = 0.99 
-6 
0.01 0.02 0.03 0.04 0.05 0.06 
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Figure 13с : Response of the LVDT for various inputs 
and directions of disolacement. 
Variation of Average Vout with Negative 
displacement. 
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уз -0.1532 + 0.2775х 8=0.99 


-0.16 
QO Column 3 
+ Column 5 
3 
0.18 Column 7 | 
(Міп з 18М) 





= - 02035 +0.5245x R=0.98 
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Рісиге 134: Response of the LVDT for various inputs 
and directions of displacement 
Variation of ( Vout / Vin) with Negative 
displacement. 
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Furthermore by comparing Figures 13b and 13d one can see that for 


the voltage region (18V-24V) the difference in slope of the corresponding 
curves is only at the third decimal point ( less than 0.5 % ). 

Consequently for this region the LVDT response is almost 
independent of the direction of displacement. 

The stability of the LVDT in time was checked by recording its 
output at a step inteval of 5 seconds for about | hour at the input voltage 
check points of 12v, 18V and 24V as shown in Table V. Plotting of the 
results in Figures 14a, 14b and 14c shows that although the output is 
pretty stable with time prior to its' stabilization overshoots for a period of 
10 to 20 seconds depending on the input voltage used (the lower the input 
voltage the longer the overshoot) and finally stabilizes rather quickly under 
a damped oscillation. It is believed that this is the reason that the data 
recorded during the initial run of the experiment showed that the 
cross-head moved up and down. Furthermore the reason that a polarity check 
did not reveal any faults, is that the amplitude of the oscillation is very 
small ( 0.0001V ) compared to the output of the LVDT which is of the order 
of several volts. 

. The check of the time dependence was made possible by monitoring 


the output of the LVDT in the neighborhood of | hour, 10 hours and 15 hours 
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TYPICAL LVDT TIME RESPONSE FOR VARIOUS INPUTS 
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Colurnn 3 
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Рідиге 140 : Timeresponse of the LVDT for Vin » 18V. 


ГІЗ 


Vout ( V ) 


-5.46 


-5.465 


-5 464 


Figure 


a 


1 


| 4c 


Vin » 24V & Pstv drctn of dsplcmnt 


G Column 4 


і 


5 10 15 20 25 30 20 40 45 50 
TIME ( sec )-Column 1 


Time response of the LVDT for Vin = 24V. 


114 


for the same displacements of 0.02, 0.04 and 0.06 inch in both directions 
and for the same input voltage check points (12V, 18V and 24 V). The results 
are presented in Tables V-A and V-B and when plotted as inFigures 15a 
and 15b showed that the minimum time dependance ( less than 1074 
relative errror ) is achieved for both directions of displacement at 24V and 
at about 0.04 in displacement. Although negligible ( less than 2х1075 
relative error ) the time dependance at 18V was almost constant with a 
shallow minimum at about 0.03-0.055 in, whereas no minimum point was 
identified in the 12V case. 

Consequently for the entire voltage region (12V -24V) the LVDT 
response shows negligible time dependence. 

Since the area of interest as far as the creep is concerned is the 
region between 0.020 in and 0.025 inches, drifting in time does not seem to 
be a serious problem for any input voltage within the above specified region; 
however a slight dependence can be observed in the lower half of this 
voltage region and hence it is prudent to use the 24V input voltage for the 


L VDTs that are going to be used. 
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TABLE V-A 


LVDT RELATIVE ERROR FOR VARIOUS INPUTS AND POSITIVE 
DIRECTION OF DISPLACEMENT 


Column 1 Column 2 Column 5 Column 4 
D 0.000030 0.001280 0.001540 


0.04 0.000810 0.001299 0.000027 
0.06 0.001320 0.001410 0.001599 


DSPLCMNT(in) OVout / Vin OVout / Vin DVout / Vin 


Dstv Drcin Vin = 12У Vin » 18V Vin » 24V 
CODES о: thr ? 10015 
TABLE V-B 


LVDT RELATIVE ERROR FOR VARIOUS INPUTS AND NEGATIVE 
DIRECTION OF DISPLACEMENT 


Column | Column 2 Column 3 Column 4 


0.06 0.00029% 0.001050 0.00 1649 
0.04 0.000519 0001113 0.000079 
0.02 0.000050 0.001460 0.001570 


OSPLCMNT(in) DVout / Vin OVout / Vin OVout / Vin 
Ngtv Drctn Vin = 12V Vin = 18V Vin = 24V 


(2025: # 415 hrs * :(Ohrs е: ihr 
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V. CONCLUSIONS- RECOMMENDATIONS 


In identifying the characterization methodology for efficient 
strength to life data generation by testing fiber bundles (each may consist 
of ten thousand (104) filaments) it was made clear that in order for one to 
arrive to meaningfull statistical data, /onmg periods of time and a very 
efficient data acquisition system is required. The high efficiency of the 
data acquisition system is primarily needed during the loading phase 
because the life prediction is based on the conditional probability that 
no failures have occurred during the loading sequence and this probability 
can be calculated only if the number of failures that occured during the 
loading sequence is Known; in the stress rupture phase, the more accurately 
the number of failures is recorded the better the results. Furthermore the 
data accuracy during the loading sequence is what gives rise to an accurate 

lower 22/7 of the Weibull life distribution which is cruisial in safe life 
prediction calculations. 

The data accuracy during the loading sequence are influenced by the 


limitations of the testing equipment and the data acquisition system. 
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The limitations of the testing equipment were examined бу а 
mechanical spring simulation and it was found that in order to avoid 
irregular indications during the loading sequence caused by the stall of the 
equipment due to internal friction, the data base of Table III-A and/or 
Figure 12c-! have to be used whereas for the 450 lbs load level of 
interest (simoultaneous testing of nine (9) graphite bundles Hercules AS4 
including the control bundle) the cross-head speed setting has to be higher 
than 0.5 in /min level with the variable speed indicator ас 21%. 

The limitations of the testing equipment during the creep-rupture 
phase were also examined by simulation and it was found that indeed the 
LVDT output is very small. More specifically, for the recommended 
excitation voltage of the 24V where the time dependence errors are almost 
inexistent a creep displacement of O. 0001 inch corresponds to a difference 
in output slightly less than 0.001V which dictates that the five digits 
resolution mode has to be used by the data acquisition system. The reason 
that the data recorded during the initial run of the experiment was 
considered unreliable (in that it showed that the cross-head moved up and 
down, whereas only upward displacement is theoretically expected from 


monotonic creep) lies on the fact that even though the LVDT output is 
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stable with time, prior to the composite sample stabilization, it 
overshoots and stabilizes rather quickly (within a few minutes) under a 
damped oscillation. This oscillation is relatively small (of the 0.0001У 
order of magnitude) compared to the LVDT output which is of the order of 
several volts. 

The data acquisition system used is not capable of recording data at a 
rate faster than 0.04 readings/sec. According to the available literature 
this result appears to be within the design constraints of the instruments 
since an interface loop involving only the HP-41CX calculator and the 
HP-3421A DACU can barely give | reading/sec. It is therefore obvious why 
the recording rate decreases substantially when the cassette drive and the 
thermal printer are also included in the loop, given that they are considered 
slow units in the first place. On the other hand the data logger routine used 
had been modified to record data at hundredths of seconds that is 50 times 
faster than the theoretical capability of the data acquisition unit itself. 
Finally due to the limitations of the HP-41CX calculator the advanced 
commands could not be used (AUTOZERO OFF, AUTORANGING, etc.) which can 
also substantially increase the recording rate. Under the above 


circumstances it is recommended to substitute the existing data acquisition 
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system with a simpler one composed of the the HP-3421A DACU and the 
HP-85 desk top computer at least during the loading sequence. The obvious 
advantages of such a system which can give a reading rate per second as 
fast as 3, 9or 11 (for 5, 4 or 3 digits of resolution repectively) are: 

l. No other peripherals that might limit the recording rates will 


have to be connected since the HP-85 has its own printer and 
cassette drive. 


16) 


More programming flexibility will be available through the use of 
the advanced commands. 


3. Substantial difference in memory capacity (319 vs16K). 


4. During the loading sequence there is no need to record the LVDT 
output due to the relative rapid changes, a fact which saves one 
channel and at the same time provides the opportunity to use the 
three digits of resolution mode (N3) which can more than triple 
the recording rate. 
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APPENDIX A 


I. INSTRON 1000 DESCRIPTION-OPERATING PRINCIPLES 


А. GENERAL 

The INSTRON model 1000 Universal Testing Instrument shown in 
Figure A-1 is set to accept a main power input of 120 +10 VAC, single 
phase, 50-60 Hz and by means of it's own power transformer, the necessary 
AC voltages are provided to the various systems of the equipment. 

All the major units of the instrument as well as their sub-units are 
described with respect to Figure A-1 by an item number given in parenthesis 


next to the unit under description. 


B. MAJOR UNITS 
These can be thought of being the loading frame (1), the drive train 
assembly and the control console (19). 
l. Loading Frame 
This is the part of the equipment where the test 


specimens are mounted and the loading (tension or compression) is applied. 
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Figure A-1: Major Units and sub units of the testing Equipment. 
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Drive Train 


This is a series of several units the prime function of 


each one can be described as follows: 


3 


Drive Motor (20) 
It's rotary motion is transfered to an intermediate drive shaft by 
means of a motor belt (21). 


Intermediate Shaft (22*) 
Its rotary motionis in turn transfered to the leadscrew drive 
pulleys (25) by means of the main drive belt (23). 


Leadscrew Drive Pulleys (25) 
These are rotated by the main drive belt (23) and are nechanically 
connected to the leadscrews (6). 


Leadscrews (6) 

These are threaded rods which are at the top supported by the 
fixed cross-head structure (7) and can turn the moving cross-head 
(2) at commanded speeds. 


Moving Cross-Head (2) 


This is the element within the loading frame (1), that applies the 
loads to the test specimen at commanded speeds. 


Control Console 


This part of the equipment serves as a housing for the 


controls and the indicator systems (19). A detailed view of it's front and 


rear panel is presented in Figure A-2. 


125 





B s — 
— = / 
5.5693 2: 














Figure A-2: Front and rear Control Panels of the testing Equipment. 
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С. OTHER IMPORTANT FEATURES 
Among the important features of the INSTRON 1000 which are 
described in full detail in the manual of the equipment [Ref. 14], the 
following two need to be generally described: 
|. Overtravel-Limiters System 
The purpose of this system is to protect the Load Transducer 
(LT, 3) from overload by stopping the moving cross-head (2). This is 
achieved in the following sequence: 
a. setting the adjustable upper and lower limit stops (13,14) in 
appropriate positions on the limit switch rod (12). 
b. Activating the overtravel limit switches (11) after the limit 


switch actuator (15) which is mounted on the moving cross-head, 
gets in touch with either one of the adjustable stops. 


D IMC Digital Processor/Controller model 700 (IMC DP/C) 


The purpose of this feature shown in Figure A-3 is associated 
with the upgrading of the INSTRON 1000 so that testing of viscoelastic 
materials could also be possible. More specifically this unit is capable to 
create a constant stress level in time, using the applied load as the control 
variable which is of prime importance during the creep. This is achieved by 
the following sequence: 
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a. setting the desired maximum and minimum load levels on the face 
of the IMC DP/C defining thus the permissible band of the load 
variation. 

D. Activating the automatic mode by selecting the appropriate button 
on the face of the IMC DP/C which allows the moving cross-head 
to move under its control to the desired maximum setting and then 
continually adjust to correct for the effects of the creep, within 
the range setting of the previous step. 

0. FUNCTIONS 
l. Load Measurement 

The load measurement is achieved by means of a calibrated 
Load Transducer ( LT, 3 ) that has to be installed into the moving 
cross-head. There are three options in calibrated LTs (10, 100 and 1000 
lbs) dependent on the anticipated load level. The input of the transducer is 
transformed to an electrical signal by means of a LINEAR VARIABLE 
DIFFERENTIAL TRANSFORMER (LVDT). This is a sensitive device that can 


measure minute deflections of the load transducer in the load weighting 


system and is composed of two major elements shown in Figure A-4: 


a. The LVDT coil assembly (2*) and 


b. The LVDT core assembly (4*). 
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LVDT coil and core assemblies. 
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It's transformed electrical signal is finally digitally displayed 
to a Liquid Crystal Display (LCD) panel located at the upper right corner of 
the front side of the console, at a frequency dependent on the commanded 
cross head speed. 

2 Extension Measurement 

The extension measurement of the moving cross-head is 
achieved by utilizing a SHAFT ENCODER (SE, 22), a device which generates 
an output in proportion to the revolutions per minute of the intermediate 
shaft ( IS, 22* ) by which it is driven. This output is an incremental 
measure of the displacement of the moving cross-head and it is also used as 
a feedback to maintain the commanded speed The defficiency of 
the extension measurement system lies on the fact that no output can be 
generated by the SE in any case case the equipment could not overcome its 
own internal friction as in a case of a partial or total stall. 

D. Creep Monitoring 

This function is achieved via the enhancement feature of the 
IMC DP/C which is capable of creating a constant stress level in time 
using the applied load as the control variable. Thus, by setting the desired 
band of load variation the moving cross head moves under the IMC DP/C 


control by continually adjusting to the effects of the creep. 
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II. INSTRON 1000 LIMITATIONS 


A OBJECTIVES 
The objectives with respect to the INSTRON 1000 limitations were: 
|. To evaluate the minimum speed that the INSTRON can be operated 
without stalling. 


2 To convert this speed to a combination of the INSTRON labeled 
controls. 


Dp To check whether this speed is load dependent or in otherwords 
given a physically verified speed check whether one is able to 
predict the corresponding stall load. 


4. To find out any other limitations that will become apparent during 
the simulation experiment. 


B. THEORETICAL MODEL 

A practical way to present a relationship between a load level (p) and 
the equipment  cross-head speed (v), is to follow the simple relations 
presented below where t represents time and І а physically measured 


displacement: 


2 


(A-1) (р/1)Ж 1/р) = 1/1 
(А-2) 1/1 = ү 


(A-3) p = v*{ 1/(1/p)]*t 


Given that the ratio (1/p) can be thought of representing the elastic 
constant of a linear model one could anticipate the following results: 
р The functions of 1 = f(t) and p = f(t) to be single valued and 
monotonically increasing. 


2 As load increases the displacement and time intervals have to 
show a tendency to decrease. 


С. MECHANICAL SET-UP 

A helical spring as the one shown in Figure A-5 was decided to be 
used for the simulation experiment provided that both the cross-head speed 
(v) and the time (t) could be easily recorded, because it was cosidered the 
best solution from the safety standpoint. 

Under the requirements of comforming with the necessary safety 
precautions for the operator and the equipment limitations, the mechanical 
set-up shown in Figures A-6 and A 7 was designed which provided the 
capability of easily achieving the target load level of the 450 Ibs order of 
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OIMENSIONS : inches 


HELICAL SPRING 


SCALE 2 





Figure A-5: Helical Spring employed for the Simulation. 
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Assembled components for the Simulation. 


Figure A-7 
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magnitude. On the other hand since one could not rely on the INSTRON 
cross-head control knob markings for the evaluation of (v) due to inherent 
limitations of the equipment [Ref. 14] and therefore a dial indicator (DI AL) 
capable to measure displacement in 0.001 of an inch ( ten times better 
resolution than that of the LCD meter ) was also used during the simulation 
procedure as shown in Figure A-8. 

The overall mechanical set-up for the simulation of the loading 


sequence is shown in Figure А-9. 


D. METHODOLOGY 


The methodology followed consisted basically of: 


|. Evaluating the quality of the response as outlined in steps B.1 and 
B.2 above. 


2. Determining the lower threshold speed that the INSTRON could be 
operated without producing irrelevant indications and the stall 
load corresponding to each one of the selected speeds, using 
the method of bracketing between selected upper and lower 
bounds. 


3. Observing the quality of agreement between the DIAL indicator 
utilized and the LCD meter. 
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Figure A-8: Dial indicator employed to measure the physical 
displacement of the INSTRON 1000 moving 
cross-head. 
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Figure А-9: Overall Mechanical Set-Up during the Simulation 
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E DETAILED PROCEDURES 

The detailed procedures used to achieve the data base available are 
presented below in a check list form to provide for ease in reproducing 
parts or the totallity of the results, under the desired incremental values 
depending on the application of interest. 

l. Familiarization 

Make sure that you are thoroughly familiar with the controls of 
the INSTRON as well as that you have already read and comprehended the 
following istructions. 
CAUTION 
DO NOT ATTEMPT TO OPERATE THE INSTRON UNLESS 


YOU ARE FEELING COMFORTABLY WITH ITS CONTROLS 
AND THE FOLLOWING INSTRUCTIONS 


2. No- load case 


a. Identify the Load Transducer on the cross-head. 
D. Calibrate the INSTRON accordingly [Ref.14 : ch. 5] 


(1) Go through steps ! and 2 of the calibration procedure 
(page 5-2 ). 


(2 If using the 100 Ib or the 10 16 capacity Load 
Transducers go through steps 3 ,4 and 5 (pages 5-2 & 
E 


(а) For the 100 ІР capacity Load Transducer go 
through step 6a (page 5-4) anyway, and NOTE 
(page 5-4) as well as steps 7b, 7c (page 5-5) 
only if necessary. 


(D) For the 10 lb capacity Load Transducer go 
through step 6b (page 5-4). 


(3) If using the 1000 Ib capacity load transducer go through 
NOTE and step 6a (page 5-3) anyway , and step 7a (page 
5-4) only if necessary. 


Set cross-head control knob to 10/.5. 


Set variable speed knob to 0.1 (first position next to MIN) and 
observe illumination of the corresponding red light. 


set physical displacement scale (DIAL) to starting reference. 
Dased on the settings of thc previous two steps make a rough 
calculation of what the maximum displacement will be according 
to the formula: 
( DISPLACEMENT ) - ( SPEED) X C TIME ) 

and set the upper limit stop. TIME has to be assumed that of the 
step 2j for the No-Load case and 3c for the load conditions. The 
lower limit stop has to be set so that it matches with the physical 


displacement starting reference. 


Set stop watch to starting reference . 
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Make sure that the appropriate record equipment (lists, pencil, 
eraser) is handy in front of you and that you have understood how 
the data recording will be performed. 


simultaneously start the watch timer and push the white UP 
button (tension) and observe it's illumination as well as, that red 
STOP button's light is turned off. 


Record physical and LCD displacements at one (1) minute 
intervals for the first ten (10) minutes, at ten (10) minute 
intervals for the first hour; if the INSTRON will stop earlier do not 
forget to record stall! time and load level. 


Push red STOP button and observe it's illumination as well as that 
the white UP buttons light is turned off. 


Push yellow RETURN button and observe it's illumination as well 
as that the red STOP buttons light is turned off. 


CAUTION 
THE RETURN SPEED 15 ALWAYS 20 in/min (500 mm/min) 
AND THEREFORE BE ALLCRT TO PUSH THE STOP BUTTON 
IN CASE ANY FAILURE OCCURS TO THE LOWER LIMIT STOP 


Set variable speed knob to 0.2 ( second position next to MIN ) and 
observe illumination of the corresponding red light. 


Repeat steps 2e up to 2] above. 


Set the variable speed knob to O.3 ( third position next to MIN ) 
observe the illumination of the corresponding red light and repeat 
step 2n above. If the INSTRON will not stal! do not exceed 2/3 of 
the Load Transducer capacity. 


Obtain plots of the physical displacement versus time with 


parameter the four different speed settings outlined in steps 2c, d 
m and o. | 
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Under load case 


Repeat steps a to i of paragraph 2. (No-Load case). 

Repeat step f of paragraph 2. 

Record DIAL and LCD displacements as well as force at one half 
(1/2) minute intervals for the first ten (10) minutes; hopefully 
the INSTRON will stall; do not forget to record stall time and 
load. 

Repeat steps 2k and 21. 


Repeat in the order given the following steps: 2m, 2e, 2f, 2g, 2h, 
2i, 3c, 2k and 2l. 


Repeat in the order given the following steps: 20, 2e, 2f, 2g, 2h, 
2i, 3c, 2k and 21. 


Obtain plots of the physical displacement and the applied load 


versus time with parameter the four differnt speed settings 
outlined in steps 2c, d, m and o. 
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|||. CREEP-DISPLACEMENT/STABILITY OF THE LVDT 


A BACKGROUND 

It is Known that for the graphite fiber bundles the order of magnitude 
of the displacement during a normal loading sequence at a constant stres 
rate R is about 2% (0.02). It is also known that the graphite fiber bundles 
are viscoelastic and therefore they creep and the total creep displacement 
is Known to be of the 0.004 order of magnitude or in otherwords 20% of the 
maximum displacement achieved at the end of the loading sequence. Due to 
the small! numbers that one has to deal with, during the creep experiment 
for a gauge length of | = 10 In it is necesary to use a very accurate device 
that will be capable of recording and/or displaying the minute differences 
in displacement during the creep. This is the point where the (LVDT) comes 
into play because the minute creep displacements of the cross-head of the 
INSTRON can be converted to voltages of the 0.1 mV order of magnitude 
which fortunately can be recorded and/or displayed using high accuracy 
multimeters. And since it is this creep-rupture behavior that can lead to 


strength-life data, the need to be familiar with the problems that may arise 
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like the one according to which no output has been recorded although the 
INSTRON cross-head had moved, is therefore obvious. Some possible 
explanations of the problem just described could be either the very small 


magnitude of the voltage output or drifting in time. 


B. OBJECTIVES 
Therefore the objective is to test the stability of the LVDT with 
respect to the following variables: 
The Input Voltage (Vin). 
The Direction of Displacement (d+) or (d-). 


The Time (t). 
The Temperature (T). 


DUN 


Based on the fact that the actual experiment will be conducted 
essentially in constant room temperature ie. within temperature 
differences of less than five (5) degrees, the stability of the LYDT with 


respect to the temperature does not appear to be significant. 


C PROCEDURE 
To implement the above mentioned objective the set-up presented in 


Figure A-10 was used which is composed of the following components: 
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l. A regulated position test stand for the mounting of both parts of 
the LVDT. The core was mounted on a micrometer's vernier by 
which a fully controlable insertion of the core at increments of 
0.001 in was possible. 


2, A power supply (НР-16168) by which the control of the input 
voltage Vin was possible. 


E: A high accuracy multimeter (FLUKE 8840A) for the exact 
measurement of the output. 

The maximum input of the SCHAEVITZ LVDT is 24V and therefore for 
time saving purposes the " bracketing technique was used between the 
limits of 24V,12V and IV respectively. If this technique would show that 
some voltage region was immaterial, bracketing would have been then 
applied for the remaining voltage region of interest. 

For positive direction of displacement the output voltage of at least 
three (3) points was obtained, for various selected values of the input 
voltage ( at least 5 ), and repetition of the same procedure for the negative 
direction of displacement was accomplished. 


The anticipated results according to the specifications were: 


|. The output of the LVDT to be independent of Vin 
Check : The ratios of Vout/ V in versus d, plot under the same 
Slope. 
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The output of the LVDT to be independent of d* or d- . 
Check : The curves of Vout’ Vin versus d+ and Vout’ Vin versus 


d- coincide. 


The linearity of the LVDT. 
Check : Curves are straight lines. 





APPENDIX B 


I. POINTS OF CONTACT 
The Hewlett-Packard engineers/technicians who served as points of 
contact in the recording rate optimization problem of the HP-3421A DACU 


and their telephone numbers are listed below for further contacts if 


necessary: 
l. Tom Turney : (408) 988-7304 
2. Russ Mc Ugh : (415) 857. 6241 


3. Ashley Werrpass : (415) 857-5673 


2. OFFICIAL CORRESPONDENCE WITH HEWLETT-PACKARD 

According to HP-iL Module OWNER'S MANUAL [Ref. 17 : p. 63] the 
following official letter was forwarded to the Customer Support Division of 
the Hewlett-Packard Co. The purpose of this letter was to explain the 
existing situation with respect to the slow data acquisition rate of the 
НР-3421А DACU, have the Company confirm the information released by 
phone about the recording rate capabilities of the equipment depending on 
the peripherals used within the HP-IL and finally have it propose viable 


solutions utilizing the already available equipment. 
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DEPARTMENT OF THE NAVY 
NAVAL, POSTGRAOUATE SCHOOL 
MONTEREY. CA 339423 5100 ім тєр ү ВЕРЕЯ ТО. 
Curriculum of Aeronautical Engineering 
Ма јог (НАР) О.М. Petridis, Code 31 
Monterey, CA 93945-5000 
November the 21st 1986 





Oear Sirs, 


| am an officer of the Greek Air Force, currently stationed at the U.S 
Naval Postgraduate School, for a Master of Science Degree in Aeronautical 
Engineering. My thesis research is related to the graphite fibers tensile 
testing, in an effort to obtain a strength-life data base that will serve in 
predicting the service life of a component made of graphite fiber 
composites, given that the strength of the fibers Is known. 


The experimental set-up is composed of a tensile testing equipment 
(INSTRON 1000) where the graphite fiber bundles are mounted. Failures of 
individual fibers within the bundles are transfered to an HP-3421A Data 
Acquisition Control Unit by means of a set of high precision load cells 
(INTERFACE SSM-100). Installed in this unit are two 10 cnannel multipiexer 
boards ( 44462A-OPTION 020 ). During the experiment, data are recorded 
using the HP-IL communications circuit as shown In the attached block 
diagram. The HP-44468A Data Acquisition Control Pac provided the basic 
data iogger software program. This program was modified to permit the 
recording of time in hundredths of seconds and edited to reduce the number 
of data registers used. A copy of this program and its output are also 
attached to this letter. Contrary to what is advertised in page 32 of the 
Operating, Programming and Configuration Manual of the 34214 the reading 
rates obtained are ої (Ле 0.04 order of magnitude (see attached plot). 


A preliminary telephone contact with your offices in Palo Alto on 
September the 29th ((415) 857-8175] provided us with the information 
according to which the reading rates reduce to NSC(I)/NA(2)/N3(2) if the 
НР-41СХ Is interfaced to the 3421A as snown in the attacned DIock diagram 
and can be substantially improved to NS(3)/N4(9)/N3(11) respectively, in 
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case the calculator, the thermal printer (HP-82162A) and the digital 
cassete drive (HP-82161A) are substituted Dy the НР-85А desk top 
computer. | 


Similar information was provided by your offices in Santa Clara 
((408) 988-7304] on October the 27th, and Palo Alto ((415) 857-8241/5673] 
on October the 3lst and suggestions to use programming utilizing the so 
called advanced commands (AUTOZERO OFF and the less digits of resolution 
possible) were also made . 


Based on the information provided so far, you are kindly requested for 
the following 


І. To confirm by letter the above information and if possible provide 
us with any insight associated with Lhe substantial reduction of the reading 
rates incase the 24214 is interfaced with either the hand-held calculator 
or the desk-top computer. 


2. To make any helpful comments on the modified HP-44468A Data 
Logger Routine and also inform us whether any faster reading rates could be 
obtained by reducing and/or by-passing a number of the HP devices in the 
loop; If this Is the case what are your suggestions for the optimum scheme? 


5. To inform us about the reasons that will improve the recording 
rate in case the HP-85A is interfaced to the 3421A via the HP-IL. 


4 What increase of recording rate can be achieved through the 
advanced commands with either one, the HP-4|CX calculator or the HP-85A 
computer. 


Your early reply will be greatly appreciated. In the mean time we 
remain, 
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PS 





Major (HAF) D.M .Petridis 


In case you have any questions please call the Naval Postgraduate 
Schoo! Aeronautical Engineering Curriculum Office at 
(408) 646-2491/2 and leave your message. 


If later than December the iSth 1986, please forward your 
answer to my thesis advisor in the following address : 


U.S Naval Postgraduate School 
Department of Aeronautical Engineering 
Dr. Edward M. Wu, Code 67Wt 

Monterey, CA 93943-5000 
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Figure A.1. Modified HP44468A Data Logger Routine (cont'd) 
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FigureA.1. Modified HP44468A Data Logger Routine (cont'd) 


37 


156 


ey НЫ а a o a 


418 
419 
420 
471 


— «= 
4 


422 
дод 
~ 
em. ы) 
426 
27 
423 
429 


RCL 25 
28 


Y 
A 
ere 15 
670 06 
"AFTER 
АКС. Х 
AVIEW 
GTO 05 
END 


Figure Al. Modified HP44468A Data Logger Routine (соп) 


38 


157 





DATA LOGGER OUTPUT FCRMAT 


O1 
02 
95 
04 
OS 
O6 
Q7 
08 
OF 
10 
11 
12 
15 
14 
15 
16 
17 
18 
12 
20 
=: 


- 
ба «ік 


ң аң» 
fn 


24 
nd 
са 


--у 


-ә 2 
23 
29 
30 
21 
-— m, 
= = 


Бэл н 
> ~ 


O 

24 
“таи 
uf ый 


29 
у 
- f 


38 
27 
40 
41 
42 
4 
42 
45 
46 
47 


LƏL DERI 
"ro 

E 

ЕБТ 08 
"I 

RCL 25 

STG O1 


LBL 10 
xm 

RCL IND O1 
ENTER 
FRC 

D ES 

+ 

INT 
Х<2Ү 
INT 

+ 

ISG O1 
GTO 10 
ЕТА О 
CLA 
АЕС. 22 
ЕСТ 09 
GTO OZ 


= AC 
г «жм 


CREATE 
Ба Е 25 
BU o 
PURGE 
CREATE 
ШЕШ 07 
"PASS " 
ДЕВЕ 22 
AVIEW 
ЕНА 
ден => 
О 

FS? 09 
SEEKR 
Ша 
DATE 
STO оо 
ЕЛХ 4 


48 
49 
50 
51 

52 
55 
«4 
== 
54 
57 
“8 
59 
50 
51 

52 
55 
54 
65 
56 
57 
68 
59 
70 
71 

72 
78 
74 
75 
76 
I7 
78 
79 
во 
81 
82 
85 
84 
85 
86 
87 
ag 
89 
90 
21 
92 
Ва 
94 


ADATE 
TIRE 
STO O1 
ЕТА 6 
ATIME 
FC? 08 
AVIEW 
FRA 

ADV 

. OO1 
ps 99 
WRTRX 

2 

STO 23 
СРЭРТО 
85-02 
RCL IND 
S10 70} 
XROM 
ASTU OO 
ASHF 
ASTO 56 
= 

зло 27 
LSL o4 
ЕТА, О 
CLA 
RET 
XEQ IND 
FS7TC 10 
ASTO OO 


eng. IND- 


CLA 

REL #01 
INT 
ARCL X 
RCL IND 
1595 27 
SIGN 


10 + “ 
з 


ENG 5 
ARCL X 


ARCL 09 
БЕЛОМ 
AVIEW 
PRA 


> 
= toed) 


“ DECODE" 


«г 
= 


ть 

97 

98 

99 
100 
101 
102 
105 
104 
105 
106 
107 
108 
109 
110 
111 
Lie 
Lig 
114 
Im 
116 
fig 
118 
119 


ГА 

1221 
-чеч 
озшш 


a= 
"Emo e) 


СРО 
ISG Ol 
GTO O4 
CLA ‘ 
FIX О 
ARCL TE 
REL 153 
mmc O9 
БЕСКЕ 
ACE IND 
cu ОД 
REE 3 
EN * 23 
1 

ee 

К 
50209 
WRTRX 
ADV 
155855 
57602785 
REL 55 
FRC 

2а 

+ 

S10 52 
RTN 
END 


2 


Figure A.2. Modified HP44468A Output Format Routine 
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APPENDIX С 


l. ESTIMATION OF POWER SUPPLIES NEEDED 
a. Load Cells (LC) Case 

The Specifications of the HP-6216B give a voltage range of 
0-25V DC and a current range of O-400mA, whereas the recommended 
excitation voltage of the INTERFACE LCs is IOV DC and their resistance 
250% 250) 

A parallel solution evidently does not work as it can be easily 
observed by Figure C-1 where the current through the loop exceeds what 
is dictated by the power supply specifications for the recommended 


excitation voltage of the LCs as shown below: 


1 = VRE RY, = v REC, 22, , (1/8) 10У/ 22(1/38.8880) 
1, = 0.514A> 0.400 A 


A series solution involving two (2) INSTRONs in which nine (9) 


LCs are connected in parallel as shown in Figure C-2 suggests that no 
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Figure C-1: Load Cell Sets connected in parallel. 
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Figure C-2: Load Cell Sets connected in series. 





more than two (2) LC sets can be powered by one (1) HP-6216B power 


supply as shown below: 


Switch ON :Z9,., ( I/R) VR. = 9 % (1/350) = 0.257 (1/0) 
із | tot 


Rtot = 38.888 0 
Баг 


Switch ОРЕ Ус › КД + V2i¢ =2* (V1 Ic) = 2 *(R tot ж ho?) 
Ус * V2i¢ = 2%38.888%0.257 V = 19.989 V 
Hence MAX Ve = 25 V>19.989 V 


However this conclusion does not take into account the possible failure of 
one LC (short-circuit). If this is the case the rest of the LCs will be 
possibly ruined and consequently it is not worthy of taking the risk. 

Therefore it is decided to use four (4) power supplies, one for 
every LC set. 

D. L VDT case 

According to the specifications of the LVDT the recommended 
excitation voltage is 24V and the highest output resistance 1090. 

Obviously in this case the series solution does not work since 
the maximum capacity of the power supply is only 25V and therefore only 


parallel connection schemes need to be examined. Under the specification 
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conditions the loop is limited by the 400mA current and hence by a total 
more than two (2) LC sets can be powered by one (1) HP-6216B power 


resistance of 60 2 as shown below: 


МАХ Io = VEGAR ot Or Rtot = Ve INAK i, = 24V/0.4A = 600 


Hence one possible yet economic situation is the one shown in 
Figure C-3 where the four (4) LYDTs are connected in parallel to the 
power supply yieling a total resistance of 60Q. According to this situation 
each LVDT amounts for a resistance of 240 Q «« 109 N and a current of 
100mA. This fact means a voltage differential of 24V which is well 
within the capabilities of the HP-6216B power supply. 

Therefore given that the LVDTs reliability is very high and 
that their wiring connections will be spatially that far appart, that the 
short circuit probability is very remote, one (1) power supply is adequate 
for all the LVDTs that will be externally mounted to the fixed cross-heads 


of the INSTRON machines. 
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Figure C-3: LVDTs connected in parallel. 
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2. IMPLEMENTATION OF THE OVERALL EXPERIMENTAL SET-UP 
The necessary procedures for the implementation of this task follow 
in a check list format. However in case clarifications are needed the 


manuals have to be consulted accordingly [Refs. 14, 15, 16 & 17]. 


a. Identification of the appropriate LTs to be installed in the moving 
cross-head of the INSTRONSs. In this case the 1000 Ibs LTs had to 
be used since the order of magnitude of the load at the end of the 
loading sequence for the test of nine (9) fiber bundles 
simoultaneously has been dctermined to be of the 450 Ibs order of 
magnitude. 


D. Removal of the existing 100 lbs LTs from the INSTRONS 
cross-heads, removal of the LVDTs from the 100 LTs and 
installation to the 1000 LTS. Installation of the 1000 LTs into the 
moving cross-heads and calibration of the equipment in accordance 
with INSTRON 1000 manual [Ref. 15], ( Figure A-4 ). 


C. installation of the Load Cells (LC)/Upper Differential Mechanism 
System (UDMS) into the moving cross-heads and of the Lower 
Differential Mechanism System (LDMS) into the lower end of the 
fixed cross-heads as shown in Figures C-4 & C-5. | 


d. Connection of the LC wirings to the specifically fabricated 
Distribution Boards (DB) by soldering ( Figure C-6a ). 


e. Connection of the DBs to the HP-3421A  DACUs and to the 
HP-6216Bs (excitation voltage of the LCs) by appropriately 
marked wiring, according to the correspodence given in Table 
МІ. 


[. Mounting of the LVDTs to the INSTRONs fixed cross-heads by 
utilizing the  specificaly designed attachement pads ІП 
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Figure C-4: Load Cells/Upper Differential Mechanism system. 
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Figure C-5: Lower Differential Mechanism system. 
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Figure C-6b : LVDTs' Distribution Board. 
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TABLE VI 


LOAD CELLS & LVOT TO HP-3421A DACU 
WIRING CONNECTIONS 


HP - 3421A DACU 


SLOT 0 (LOWER BUS) SLOT 1 (UPPER BUS) 
CHANNEL * CONNECTED ITEM CHANNEL * CONNECTED ITEM 

0 Actuator 0 Load Cell*8 

| Actuator | Load Се 19 

2 EC Exctn Vitge 2 LVDT 

б Load Cell*| j LVDT Exctn Vitge 

4 Load Cell*2 | 

5 Load Cell *3 

6 Load Cell*4 

7 Load Cel|*5 

9 Load Cell*6 

9 Load Се!1=7 
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predetermined locations so that no negative voltage outputs could 
be obtained that could delay the data acquisition rate ( Figure 
C-7 ). 


Connection of the LVDTs wiring to another specifically 
fabricated distribution board (db) by soldering ( Figure C-6b ). 


Connection of the db to the HP-6216B (excitation voltage of the 
LVDTs) and to the HP-3421A DACUs by appropriately marked 
wiring according to the correspondence given in Table МІ. 


Connection of the HP-85 desk top computer to the HP-3421A 


DACUs through the HP-IL for the monitoring of the loading 
Sequence. 
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LVDT to INSTRON 1000 attachment pads. 


Figure C-7 
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